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ABSTRACT 


In this work the coaxial waveguide antenna is treated by the 
Wiener-Hopf technique and the transient radiation from a cylindrical 
monopole is developed in the light of the rigorous results obtained 
from the Wiener-Hopf analysis. Analytic expressions are derived for 
(1) the electromagnetic fields in the feed line and (2) the far zone 
radiation fields of the coaxial waveguide antenna, with time harmonic 
excitation voltage. Complete characterization of the transient beha- 
vior is also found for (1) the fields interior to the feed line and 
(2) the radiated fields for excitation voltages arbitrary in their 
time dependence to the extent that kb,ka << I . This corresponds 
to the case of a thin antenna and excitation voltage with a non-zero 
rise time, specifically chosen so that frequencies violating the 
restriction kb,ka «« ] are negligible. 

The transient radiation from the cylindrical monopole is 
developed in a closed analytic form which is relatively easy to 
interpret and apply. The expressions found offer an alternative to 
transient analysis by conventional methods requiring numerical tech- 
niques involving extensive computer calculations. They are also the 
basis for an uncomplicated procedure to synthesize a desired behavior 


of the transient radiation from cylindrical monopole antennas. 
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1. INTRODUCTION 


This dissertation presents an analysis of the coaxial waveguide 
antenna by the Wiener-Hopf technique and the determination of transient 
radiation from the cylindrical monopole antenna. This work is the 
first successful treatment of a semi-infinite coaxial structure by the 
method of Wiener-Hopf and the first rigorous determination of transient 
radiation emanating from an open ended semi-infinite waveguiding struc- 
ture. In addition, a new and simple method for calculating the 
transient electromagnetic radiation of a cylindrical monopole antenna 
is presented. 

The study of electromagnetic boundary-value problems by the 
Wiener-Hopf technique began in the early 1940's. Among the most notable 
of the early workers were Schwinger [1] and Copson [2]. The early 
applications of this powerful method were inthe rigorous solution for 
the fields within a bifurcated parallel plate waveguide and diffraction 
of plane waves by a conducting half-plane. The method was later used 
to solve the problem of radiation impinging on or emerging from the 
open end of an unflanged circular pipe and a pair of parallel plates. 
Much of the early history of application of the Weiner-Hopf method is 
detailed by Bouwkamp [3] and an authoritative account of the method was 
published in 1958 by Noble [4]. In this work the integral equation 
approach of the early workers is not stressed. A more direct simple 
method due to Jones [5] is emphasized. It is the application of 
Wiener-Hopf technique due to Jones that will be employed within this 


report. 





d 
There have been several recent investigations of the transient 


behavior of antennas [6-15]. The focus of many of these investigations 
is on the time behavior of the radiated electromagnetic field. Refer- 
ences [6,8] and [14] present theoretical as well as experimental 
results. In most of the prior research efforts, the problem is formu- 
lated in terms of circuit theory and the antenna is treated as a lumped 
circuit element. The feedline and the exciting sources are replaced by 
their Thevenin or Norton equivalents. The formulation of the problem 
in terms of circuit theory yields complex frequency domain fields which 
must be Fourier transformed if the time behavior is of interest. Such 
calculations require numerical methods due to their extreme complexity. 
The work of Schmitt, Harrison and Williams [6] is a determina- 
tion of the transient radiation and reception performance of a thin, 
finite length cylindrical monopole over a perfectly conducting infinite 
ground plane with a coaxial feed. In the circuit analysis formulation 
(frequency domain) the antenna was represented by an impedance and the 
radiated field related to the voltage across the antenna impedance 
through a proportionality factor, called effective height. The time 
dependent behavior is obtained by Fourier transforming the frequency 
domain field. A similar procedure was used to obtain the properties of 
the antenna in reception; the induced voltage on the receiving antenna 
was related to the incident electric field intensity through the effec- 
tive height. The theoretical work of Schmitt, Harrison and Williams 
was supplemented in 1971 by Abo-Zena and Beam [9] to include the time 
behavior of near zone field intensities and in 1970 by Palciauskas and 


Beam [10] to include the far zone time behavior at latitudinal angles 





de 
other than 90°, which was the only angle considered by Schmitt et al. 


King and Harrison [8] formulated a frequency averaged reflec- 
tion coefficient to theoretically obtain the reflected electromagnetic 
fields on the coaxial feed line of an infinitely long cylindrical 
antenna or for the time interval before the excitation pulse can 
reach the tip of a finite length antenna. Use of frequency averaging 
is necessitated by the fact that the antenna input admittance is 
frequency dependent, which constrains the conventional transmission 


line reflection coefficient to be frequency dependent. 


The transient current on a circular tubular infinitely long 
antenna that is excited by a voltage which is a step function of 
time was determined in 1961 by T. Wu. He later (1969) furnished a 
corrected analysis for this current in Reference [11]. An earlier 
work (1960) by Brundell [12] treated the more general problem of 
determining the space and time behavior of the electromagnetic field 
for such an antenna. This work was also supplemented in 1970 by 
Latham and Lee [13] who found the early and late time behavior of 
the field radiation from a hollow infinite cylindrical antenna 
excited by a voltage that is a step function of time. 

Lamensdorf [14] determined that the time behavior of the 
coaxial cone antenna in reception is primarily a time derivative of 
the incident electric field intensity falling on the aperture. This 
work also contains a heuristic discussion of the time dependent elec- 
tric field radiated by a cone antenna. Chang [15] analyzed the tran- 


sient reception characteristics of the annular slot antenna using the 





M. 
circuit theory approach. References [14] and [15] are recent works, 


1970 and 1971 respectively, and have application as part of a receiv- 
ing system for measuring scattered impulse responses of various 
objects or as a receiving antenna on an aircraft or guided missile. 

Inthe systematic study of radiation from the coaxial aperture 
antenna which follows, analytic expressions for the electromagnetic 
fields of the coaxial aperture antenna are given. These expressions 
account for waves traveling in both the forward and reverse directions 
on the feedline and the radiation fields for harmonic time dependent 
excitation. Complete characterization of the transient behavior is 
found for the fields interior to the feedline and the radiated fields, 
for excitation voltages which are arbitrary to the extent that 
kb,ka << 1 . This corresponds to the case of a thin antenna and 
excitation voltages with non-zero rise times, specifically chosen so 
that frequencies which violate the restriction kb,ka «« ] are neg- 
ligible. 

The annular slot antenna is considered to be a special case of 
a cylindrical monopole antenna (h=0). The fields radiated from the 
annular slot antenna are approximated by well known methods in which 
the incident field on the aperture is considered to be the aperture 
field. The radiated fields so determined are interpreted in light of 
the rigorous results obtained from the coaxial waveguide antenna and 
a model is obtained for a cylindrical monopole antenna in which the 
height of the radiating element is greater than zero (h > 0). The 
transient radiation field intensities are obtained from this antenna 


model. The determination of the transient radiation through the 
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procedures developed in this research offer an alternative to comput- 


ing transient behavior by numerical methods where extensive machine 
calculations are required. The procedures reported on herein can also 
be used as the basis of an uncomplicated synthesis procedure to produce 
a desired time behavior in fields radiated from the coaxial waveguide 
and the cylindrical monopole antennas. 

The results of this research should prove to be very useful in 
the areas of 1) determination of antenna properties (driving point 
impedance and radiation pattern) through time domain measurements [16], 
2) identification of and discrimination between conducting bodies 
through radiation scattering from their surfaces and 3) the calibration 
of equipment designed to measure the effectiveness of shields designed 
to protect equipment or installations from intense, short rise time 
electromagnetic pulses. In each of the above areas a precise knowledge 
of the time-behavior of the radiated electromagnetic field is required. 
Additionally the ability to synthesize a particular time behavior in 


the radiated field is highly desirable and very important. 
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2. THE COAXIAL WAVEGUIDE: AN APPLICATION OF THE WIENER-HOPF TECHNIQUE 


In this chapter the frequency (spatial) domain electromagnetic in- 
tensities for an open ended coaxial waveguide will be derived. This will 
be done using the methods of Wiener and Hopf [4]. Use of this technique 
will provide a determination of transverse electromagnetic mode (TEM) 
as well as higher order mode intensities within the structure. 

We take the coaxial waveguide structure to be that of Figure 2-1. 
The inner conductor is a hollow tube. Both outer and inner conductors 
are assumed to be semi-infinite in length, infinitesimally thin and 
perfectly conducting. The cross section is circular providing for a 
coaxial annular region between the ccnductors. We also assume that the 
waveguide is embedded in a homogeneous, isotropic, dielectric medium 


having the constitutive parameters of free space. 


Specification of the Problem 


We begin the specification of the problem by listing the Maxwell 


equations: 
ZA 
VxE = ” 
> -> D 
VxH = Jt E 
Uv. D - p 
er (2.1) 


The fields which satisfy the above equations must also satisfy the con- 


stitutive relations: 








Figure 2-1. The coaxial waveguide antenna 
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Under the assumption that the coaxial waveguide has been excited in 
TEM mode, there is a traveling electromagnetic disturbance proceeding 
in the positive z direction. The time origin is taken to be the 
point in time when the leading edge of the TEM mode exciting field is 
incident upon the aperture located in the z= 0 plane. 

The specific method used to initially excite the coaxial struc- 
ture and the field distribution prior to the time origin are not of 
interest in this work. 

It is well known that the TEM mode field of a coaxial waveguide 
is independent of angular variation [17, p.326] in the cylindrical 
coordinate system. The coaxial structure also does not introduce any 
boundary conditions which have angular dependence. From these observa- 
tions we conclude that the scattered field must also be characterized 
by the same symmetry (i.e., independent of angular variation in cylin- 
drical coordinate system). 


In cylindrical coordinates the field components can be written 


as 
É PE a, + Ela 
= a a a 
0 p > > Z Z 

and 
> -> > > 
H = Ha_+H,a, + H_a 


-J 
where a ,a 





-98 
We now Fourier transform the Maxwell equations to change from a 
time domain to a frequency domain representation. The Fourier trans- 


form pair [18] is defined as follows: 


ει | f(t) ett at 
ver } 
f(t) = — | N Sm 
ver 2 
When the first two Maxwell equations are multiplied by L gut and 


Vem 
operated on by the integral operator, the results are: 


Lor tete dt=-—L EDT 
ven} ven ; 
and 
ος = E i 
A | v S lot dt - C (+ elot dt 
ven 2 von ὁ 


Interchanging the curl operator and the integral sign on the left hand 
Side and integrating by parts on the right hand side yields equations 
in the frequency domain. To form these results, use has been made of 


the fact that the field intensities must be zero as |ζ| » ο. 


Vx S = + iw Bla) 
V x H(w) = J(w) - iw D(w) 


It similarly follows that 





ES 


Ve D(w) 


" 
m 
— 
e 
indi 


ν «Βίω) - 0 


B(w) = uw) and D(w) = o E) 


When the field components in cylindrical ooordinates are substituted into 
the frequency domain Maxwell equations for a source free region (i.e., 

> 

J(w) is set to zero), we obtain the following relations between the 


field components: 


ΘΕ 

iw B. (2.1.1) 
ΘΕ JE 

iw B, = + (22) 
JE dE, 

Iw By = TAIN (2.1.3) 

and 

= sone 

-Ίω ae - az Ma (2.1.4) 
H 9Η 

-1W D, = et E (2.1.5) 
oH ΘΗ, 

-10 P = Dm (2.1.6) 


By using the constitutive relations and combining equations (2.1.3), 


(2.1.4) and @.1.5),we derive the scalar wave equation for H 
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This equation, upon simplification, becomes: 





9 12 1 9 
[ + = — - — + —+ ] Η (ω) = 0 (272) 
ma D ορ o μα. E $ 
where ue, =t e ean T Bl a e 
€ 
produce 
2 2 2 
9 la 1 9 ω 
[— *——-—*—t-—xJEQ(w = 0 ο 
oot p IP ot 37° E $ 


From the foregoing equations it is apparent that the solution of 


(2.2) and (2.3) for Hy and Es is enough to completely determine all field 


components. The work incorporated in Appendix A shows that Es mus t 


be zero for the coaxial waveguide, therefore the only field intensity 


components present are Ha E, and E, . The rest are identically 
zero. The analysis has been reduced to solving the boundary value 


problem 
2 2 
9 9 ep 2 
[— + *t-—-—tk]lH,(o,z,) - 0 (2.4) 
m "a o 3p 2 $ 
H (0.25) nis -0 < Z <% (2.4.1) 
E,(a,2,u) - E, (b. z,0) = 0; z<0 (250723 


The Sommerfeld radiation condition, 


9Η 
lim (> - ik H,) = 0 must also be satisfied. 


r> œ 


$ 
At this point it is useful to point out that in the above equations 


2 ut 
k = = (2.4.3) 





Eno 


is a scalar quantity and equation (2.1.5) can be written as 


-iweg, = Lá [o n (2.5) 


This will be of use later. 

The well known expressions for the TEM mode fields for the 
coaxial wavequide [17] will be given here to set the notation. We 
take the field intensities incident on the aperture from the negative 


z direction to be 











pi _ _V(w eikz (2.6) 
P gnbya P l 

; E ikz ikz 

M [ο víw) e x e 
un = m in b/à er M(w) $ (2.7) 


v(w) is the frequency domain voltage between the inner and outer con- 
ductors. 

The space containing the waveguide is divided into three regions, 
portrayed in Figure 2-2. Equation (2.4) applies to each of these regions. 
The field satisfying this equation in region one is designated 
Hg (02.0) . Similarly we have Hs, (paz) for region two and 
H gi 9290) for region three. The equation for each of the regions is 
a partial differential equation in the variables p and zZ. By 
Fourier transforming in the spatial dimension zZ , an ordinary differen- 
tial equation in the variable p is obtained. 

Following the notation of [4] , the spatial Fourier transform 


pair is defined as follows: 
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Region | 


ΜΑ. ae LL Region ZZ 
P^ Region 3 A 





2b 


BLEI. 


Region 1 


[A 


b<p 
Region 2 a p NW 
p<a 


Ja JA 


Region 3 


Figure 2-2. Division of the space surrounding the coaxial wavegui de 
antenna into regions 
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Fla) = — | f(z) οἱ dz 
νοπ 200 

flz) = E | F(a) e 12 da 
von «200 


The "plus" and "minus" functions are defined as 


where 


Fla) = Lf f(z) e'™ dz 
0 


Jon 
0 
F (a) = ΡΝ 
var «OO 
F(a) = F(a) + F(a) 


εις ττ is complex. 


(2.8) 


ο δη) 


(28:2) 


(258:3) 


(2.8.4) 


Using well known theorems involving Fourier transforms in the 


complex plane [4], the regions of analyticity of F(a) and F (a) 


will be found. 


Suppose that we are given 


numbers. In addition suppose that 


[f(z)|] < A exp(t_z) 


where A is a constant. Then 


|F'(a)| « 


where M< œ 


f(z) e!?^ qz 


an 








1 
Ven |, 


as 


Es 


von 


EE S, all real 


is chosen large enough for the inequality to be valid. 
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The last integral of the above expression converges if and 


only if 


IZ 


or 


By similar reasoning, it can be shown that if |f(z)| < B exp(t,z) 


as z > -œ (B constant) then ΙΕ (a)| is convergent provided 


< 
T Cy 


Therefore pita and F (a) are convergent in the strip 
T_<T<T, Of the a plane. 

We also note that if IF" (a)l and |F (a)| are convergent 
in a region of the complex plane, then the magnitude of nth ordered 
derivative (n < œ) of a) and F (a) with respect to a also 


exists in that region. In fact such functions are nth order continu- 


ously differentiable in the region. 


d IF (a)! e 109) 














d ‚+ 
στι (F in E 





it 
x τ. 


< + ze dz 


M 
. Il > 
i] z"f(z) e!97 dz 
Vm y 








| 
Ver y 





Er 
where M«* is properly chosen to ensure the inequality.  Concen- 


trating on the term of significance, we write 


E (ia*1 )z (ia *c )z[ » (ia* t )z 
| z" e " dz= Ze s - ZI e τ 
M M ees 





dz 


After n integrations by parts, the above becomes 


Y n! UE e 12 = : y aves {oe co 
pe] ^ (iat IF M (ia +r)” M 


which is finite for finite n provided 


TEST 


Again by a similar discussion, it can be shown that F (a) is nth 
ordered continuously differentiable in the open region rt « t, provided 
να) e Bets ας 2» -ο. 

Since all that is required for analyticity in a region of the 
complex plane is that the first derivative of the function exist in that 


region, we are assured that F,(a) and F (a) are analytic in the 


strip τ «τ-«τ, Of the a plane provided of course that 


eZ 
[f(z)| <Ae Z >o% 
and 1,2 
If(z)] «Be z >-© 
Upon multiplication of equation 2.4 by l elaz and operat- 


von 
ing with the integral | dz , the following equation in 


oo 


p emerges. 





>17. 

"ο . 2 
- | E, Ha (0, 2,u)e ^" t [5 + u 
ven + dz Ip E 


dr 


1 2 B 
30 = E K J Ha (95050) zu 


Integrati ng by parts twice and imposing the conditions that HCP Zw) 
and SH (p,z,w) >0 as jz| >œ% , produces an ordinary differential 


$ 


equation in o , 


ld 1 2 
Er "ES Y] Hato 40,0) = 0 (2.8.5) 
p p 
ce (2.8.6) 


It is noted at this point that the incident magnetic intensity does 
not satisfy the differential equation 2.8.5)since it does not converge 
in the limit as z *-. This will be explored in detail a little 


later. 


For region l|: 


2 
d l d ] 2 
[ ne nn Y ] H (0,0,0) = 0 (2.9) 
do? 04» ϱῦ $1 
Hy (esau) + 0 as pco (2.9.1) 
d Lo Hy (Peau) J = 0 (2.9.2) 


p=b 
Equation (2.9.1) follows from the Sommerfeld radiation condition and 


equation (2.9.2) follows from equation (2.5) and the boundary condition 


(2.4.2). 





Zee 


For region 3: 


2 


d ld 1 2 

—>+t---—-y]H (pw = 0 (2.10) 
do” e do of $4 

Hy ο ανω) = 0 (2.10.1) 


tl 
O 


t Lo Hg (070,01 (Zee) 


p=a 
Equation (2.10.1) and (2.10.2) result from imposing (2.4.1) and (2.4.2). 


For region 2: 
In this region the magnetic intensity will be written as a 
linear combination of the incident TEM magnetic intensity and another 
term vy, (p,2,0) . 
$2 
a Z, ΖΗ 
The term Y, (0,2,0) represents obviously all of the reflected fields that 
2 
exist within the feedline, on z«0O. For the other extreme of 


2.205 Y, (0,Z,w) must account for the fields which exist in free space 
2 


from sources within the feedline. It must be representable as 


Yo (p,2,0) 7 (A + other field terms) , z>0 
2 


since Ha cannot be present in the half space z> 0 because it vio- 


lates required boundary conditions. 
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Since the intensity Ho (p,z,w) must be a solution to the dif- 
2 
ferential equation (2.4) the terms Ho and Y, (0,2,0) must satis fy 
2 


the following: 


2 2 


9 9 19 1 24 111 
[— + E +t-—-3+kJ]H = 0 
azl 30° ϱ ὃρ "d $ 
and 
2 2 
9 9 19 ] 2 
[Ez + E + — nN at k] ψ, (zw) = 0 
z^ 30° o Ip pe $5 
The incident field Ha = Ho ekz identically satisfies the 


first of the above equations. Therefore we only need to solve the 
second equation. The magnetic field intensity for region 2 will be 
completely known after this is done. 


The above differential equation for Vg, (0 22 ui will now be 


jaz 


transformed by multiplying through by Ee and operating wi th 


vor 
| dz 


2 
9 la 1 2 

πο προ ο ο ες ρα aad 
30° 9 3e p $5 


iaz 


v |% 
N m 


| 


Integrating by parts twice and requiring that UR (0,Z,w) and 
2 


ο 


AN 


dz s, (ore) +0 as |z| » , we derive the differential equation 


for region 2. 


(> + 12- L - Y] Va ίραω) = 0 (2m2) 
dp d D 2 
d 
—(o V, (o,a,0)) = 0 ου 
do" ὀρ ο 
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The condition expressed in (2.12.1) follows when it is recognized that 


d a= 
mi Ha) 0 (2.12.2) 


and use is made of the boundary condition expressed in (2.4.2)and (2.11). 
It was previously stated that the incident magnetic intensity 

does not satisfy equation(2.8.5). In the derivation of (2.8.5) it is 

essential that the function and its derivative vanish as |z| > œ 

It is obvious for slightly lossy media that Η does not vanish as 


Z > -© „ therefore it cannot satisfy equation (2.8.5). 


Upon multiplication of equation (2.11) by LL ¿102 and operat- 
2T 
ing with | dz , the following results: 
1 OZ 
H, (p,0,0) = —— | (H,) e dz *y, (0,0,0) 12123) 
ὀρ Ver ? $2 


Substituting the value of H' in the above integral produces: 





$ 
5 ikz . 
Mio) | e ϱἷαΖ gz (2.12.4) 
ν2π p 


This integral does not converge. If it is assumed that the 
medium is slightly lossy, we can write the quantity k as a complex 


number 


k = K1 + ik, ; κ] >> Κρ 5 SELS PRU) 


The loss term Ko will be allowed to go to zero later in the problem 
to yield results for a lossless medium. The value k is specifically 


chosen in this form to ensure convergence of the TEM magnetic intensity 





E 
as z> +o. This choice of k ensures convergence of(2.11.4)at the 
upper limit (z ^ +~) . Therefore it is concluded that we should operate 


on (2.11)with the integral E rather than | dz . When this is 


done, the following equation, which will be used later, results : 


iM(w) 


Hy (0,0,0) Et 
2 vet pla +k) 


+ y, (p,a,w) (2.12.5) 
$5 
In the differential equations for the various regions of the 


problem we have made use of a quantity 


The square root is 


y = Joo - k ο 


which is a multivalued function of the complex variable a defined on 
a two sheeted Riemann surface. The value of y must be uniquely 
defined so that only the proper branch of this double valued function 
is used. The branch cuts for this function are determined in Appendix 
B and depicted in Figure 2-3. 

The necessity of satisfying the Sommerfeld radiation condition 
requires that the asymptotic frequency domain behavior of the magnetic 
intensity for the space around the coaxial waveguide antenna (exclusive 


of regions 2 and 3 with z < 0) be that of an outward traveling or 
evanescent spherical wave. 
ik|r] 
Hives 


"Jn 
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ΠΠ] ΤΣ 


If the value of z is such that z >> p the asymptotic behavior 


specializes to 








| | otklzl 
H Pala : Z ^ to 
$1 |z| 

ekz 
υ ο σω © Ζ > +o 
9 Z 

elKZ 
Hp (02220 = ; Z > +o 


In region 2 with z « 0 , the magnetic intensity is known [17, 
p.327] to consist of the TEM mode and higher order E modes. Only 


Eon modes are present due to the symmetry requiring no variation with 


the angle $ . The only field intensity components allowed are Ho : 
Es and E, : 


From equation (2.11) we have 


g Carikz 
Ya Pazo) dE 
izek - (yo) 


Con {4 (Von? orsa) - te 


Ma 8 


n 
whe re ΝΙΝ (γ' É = i/(y' γέ. Ke ns lea 
Yon Yon , 9C » J5 


Y is the nth ordered root of the auxiliary equation 


Jg Y oO b? Y (Yon) = VOY D) J (Yon?) 
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C and Con are constants relative to the spatial variables 


o and z. Jy(z) and Y,(z) are the ordinary Bessel function of the 
first and second kind, respectively of the vth order. 
In region 3 with z < 0 , the magnetic intensity is also known 


[17, p.322] to be the En modes of a cylindrical waveguide. 





u 2 
oo P Uon iz κέ. --η 
Hg. pozo) m dà er At a ) e 
"Uo Mns. . _ 
where ./ k- - (e = q FE EU : ο ου”... ee C is a 


constant relative to p and z and us is the nth ordered root of 


the auxiliary equation, 
Jo Qu) = 0 


Given k = k_+ ik and the asymptotic behavior of the field in 


] Ὅν 
region 1, it is apparent from the earlier discussions that He (p,0,0) 
1 


is analytic in the strip 


-ks «ως < Κρ 


of the complex a plane. 
For regions 2 and 3 we must consider two cases, the first case is 
when the Eon modes within the region are evanescent and the second 


case is when they are not. 


Case |]: k < Yon region 2 
u 


k < = region 3 
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Case 2: k > Yon region 2 
“on 
k > πα region 3 


For Case ] the E modes in regions 2 and 3 are evanescent, 
therefore the asymptotic value of the magnetic intensity as z > -œ 


is the TEM mode portion of Ya within region 2 and zero within region 
2 


3. We conclude that in the event of Case |, 


Y, (Psu) and Aj ero) 


are analytic with the strip -k> ES Κρ of the a plane. 
For Case 2, the Eo modes are not cut off and the asymptotic 


forms behave as 


g^ 8 
-izk /1 - =y -iz(k,* ik) M $ 
e S e k , Z> -œ 


This behavior produces a strip of analyticity that isat least qiven by 


-k <T<k, Re(,/1 - (8) 


The symbol Re means the real part of the quantity that it precedes. 
Although at this point there is no apparent difficulty with this result, 


the inclusion of values of ka and kb such that 


k = u 
a on 
and 


kb = Uon 
Be 1 : 
where ME the mth ordered root of J (u) J fe u) =0 requires 


special attention. It will be shown later, that the field intensity 
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outside the waveguiding structure is zero at these frequencies. 
In light of the above, we restrict ourselves to the fundamental 


mode by imposing the conditions 
kb << ] and ka << ] 


when we are considering a time variation other than harmonic time 
dependence. Under this restriction Case 2 does not occur, and 


the functions 
Πρ (ρνανω) ; Moa and Hg, Vos] 


are analytic in the strip 
-k> ers κ, 


of the complex a plane, see Fig. 2-4. 
When we consider harmonic time dependence, the strip of analy- 


ticity is changed to 


"Κο «τε Κρ πι - (Ey 


where 8 = Yi or ot for region 2 or 3, respectively, n21,2,3,---. 
The value of B is the magnitude of the cutoff wave vector of the 
highest order ae mode which is non-evanescent in region 2 or 3, as 
appropriate. Also ka # Uon and kb # Um’ 

To obtain a unique solution to the problem, it is necessary to 
determine the algebraic behavior (as a+) of the factors in the 
completely factored and decomposed Wiener-Hopf equation [4, p.37]. 


This determination is made by considering the edge conditions [19]. 





sa) = pw 





a plane 


¡ko 


| 
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For the infinitesimally thin conductors utilized in this problem, the 


edge conditions are 


E + 
E, (0,2,u) 5T 1/2 as Z at p=aandb 


Equation (2.5) states that 


E (p,Z,w) ac 5 35 Lo Halo 2,0) ] 


The point of interest is that 


1/2 


= [o H4(o,2,u)]« z' (2.14) 


+ ; 
Bst at o = aorb. 25 0° means that z= 0 is approached 


from the direction of positive values of z. 


Development of the Wiener-Hopf Equations 


In this section the Wiener-Hopf equations will be developed. The 
solution to these equations will lead to the spatial transform of the 
desired field intensities. 

We begin by writing solutions to equations (2.9), (2.12) and 
(2.10) in that order. 


Ha (p,a,w) = Α(α,ω) Κι(γρ) (2.15) 
] 

Vo (o,a,w) = Bla,w) 1,(yo) + Claw) K- (Yp) (2.16) 
2 

H, (0,00) = D(a,w) 1, (ye) (217) 
P3 


The form of the solution to (2.9) expressed in (2.15) was chosen so 


that the asymptotic behavior would satisfy (2.9.1), thereby satisfying 
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the Sommerfeld radiation condition. 

The solution to (2.10) expressed in (2.17) does not include the 
first order modified Bessel function of the second kind since it 
becomes infinite as o > 0 a condition that violates equation (2.10.1). 
The function I, (yo) +0 as p +0, 50 ΛΝ = 0 and (2.10.1) 
is satisfied. 

The coefficients A(a,w), B(a,w), C(a,uw) and D(a,w) are un- 
known at this point. It is the task of finding these coefficients 
which leads to the Wiener-Hopf equations. Each of the coefficients 
will be expressed in terms of two heretofore unmentioned functions, 

η (α.ω) and H* (o0) , which will be found by employing the Wiener- 
Hopf procedure. To this end we will make use of the relationships 

which exist among the magnetic intensities and their derivatives (with 
respect to p ) evaluated on the boundaries separating the three regions. 


At pop = a we have 


H, (a,z,w) = H, (a,z,w) for 20.00 
$2 $3 
Therefore, 
H? (aaow) = H* (a,o,0) (2.18) 
Po $3 


From the analogous condition on the magnetic intensity at p =b for 


z > 0 , we have 


+ 


mo (2.19) 


H* (b,a,u) = H 
9 


These equations relate the magnetic intensities at the bound- 


aries and will be used at a later stage in the development which leads 
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to the Wiener-Hopf equations. 

It was assumed earlier that the coaxial waveguide antenna struc- 
ture was made of a perfectly conducting material. Such material does 


not support a tangential electric field, therefore 


E_(a,Z,u) = 0 
and 


ll 
e 


E, (b. z,u) if z< 0 


E (0, 2,u) is a continuous function of p , thus E (0,z,u) has the 
same value at a given value of o regardless of the direction of 
approach, whether from smaller or larger values of p. 


Equation (2.5) states that 


la 
ott 
ES Lm (o p) 


When this is combined with the preceding discussion about the tangen- 


tia] electric field, the results are 
$ [oH, (o,z,U)] = Š- [oH, (p,2,0)] 
dp 94 do $5 


at po =a. By a spatial Fourier transform (see equation (2.8)) of 


these equations, it is found that 


d d 
dp Lon, (mau) do let, (am). 


u 


When we consider that E, (a,z,u) 0 for z «0 , this equation re- 


duces to 





xo 


T [οὔ] (ονανω); (2.20) 





- ἃς [ού] (0200) 





p-a p=a 


Note that we have used equations (2.11) and (2.12.2) to obtain (2.20). 


By very similar arguments with p=b we also find that 


= Coll, (0.040) CN 





= a Low, losa) 





ϱΞΡ p=b 
We now define W” (a, 0) and "MCN as follows 


W law) à S 





eo) E d; Lot, loca) 





These functions of a are analytic in the upper half plane given by 
To “ko ; 

Returning to the spatial frequency domain field solutions given 
by (2.15), (2.16) and (2.17), we multiply each of the equations by p 


and take the p derivative of the resulting equations [20, p.376]. 


d i τα Ἱ-- 
d; Lem T [oH] yp Alaw) K (yo) (2.2220) 


x Low, 1 + T Low, 1 = yplB(a,w) I (yo) - C(a,u) K (vo)] (2.22.1) 


Pot + © fon” ) = 2.22.2 
D Loy ET Lu } γρ D(a,w) I (ye) ( 
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By taking the last of the above equations and evaluating it at p=a, 


it is seen that 


d 
do CN - W'(o,u) - va D(o,u) I (va) 








p-d 
or 
μ᾿ ανω 
D(a,w) = - (2.23) 
ya I ya 
Note that 
Be H. = 
do Lo al o=a : 
since, 


Recall that 
ΠΗ, (ρ.Ζ.ω)] - -- [ρη, (ρ,Ζ»ω)] - -- [ρψ, (ρ.ζνω)] 
dp $3 > 3 dp $5 3 3 dp $5 > 3 


at p=a. (This may be easily seen to stem from equation 
(2.12.2) and the continuity of the tangential electric field E). 


Since the tangential electric field is zero at p-aorb for z<0, 


| 
© 
-h 
O 
-- 
N 
N 
© 


d 
Lov, (ρ»Ζ»ω)} 
do ^", 





p=a,b 
and 


| 
© 


d = 
Wr [oyp (P00) ] 





p=a,b 
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Combining this result with equation (2.20) and equation (2.22.1) 


evaluated at p = a leads to the following: 


W'(a,w) = ya[B(a,u) I (ra) - C(a,w) K (ya) ] 


or 
W'(a,w) - yaB(a,w) I (va) 


ge Ya Kc (Ya) 


(2.24) 


Substitution of the values of D(a,w) and C(a,w) given by equations 


(2.23) and (2.24) into equations (2.17) and (2.16), respectively, yields 


+ 
_ M (o,0 


W" (0,0) - B(a,w)ya I (va) 


ὑρ,(ρνανω) = B(0,w) 1,(yp) - Ray a a Κι (γρ) 


and 


Evaluation of the immediately preceding equations at p = a and use of 
the Wronskian relationship produces the following intermediate expres- 
sions: 


Ν i W (asw) I4(ya) 
Πρ (ανανω) als ο ο. = γα Τα] — 


W (0,0) Ky (va) 
+ s . δω ; ] 
Y, la saul + Hg, asco) "va Kola] " va Ray 


By subtracting the lower of the two equations immediately above from 


the top equation, we find that 
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Hs (a,a,w) - V. (a,a,w) + H, (a,a,w) - y, (a,a,w) = 
$3 > > $^ > > $4 3 3 Yo, 3 > 


W (aw) B(a,w) 


(ya) I (ya) K (Ya) ο 


The Wronskian was used to achieve some simplification in the above 
equation . 


If we evaluate (2.12.5) at p= a we find 


u* (a,a,w) = Vy. (a,a,w) = . iM(u) — 

$5 $2 Von a(otk) 
But it is known from (2.18) that 

+ + 

9 9 = H 9 9 

Hg, (a ανω) $^ α«ω) 
Therefore, 

HT (a,0,0) πι V. (a,a,w) = — Mo) _ 

P3 $2 VE a(o*k) 
and it follows that 

Tm 
iM(w) W (a,w z B(a,w) 


+H, (a,a,w) - Y (asaw) = 
2 


/2v a(atk) P3 (ya) I (ya)K (ya) ya K (va) 


(2.25) 


Evaluating (2.22.0) at p = b we see that 


Mas) = - Fo (2.26) 





τας. 
Note that the boundary condition E „bs Z "S. = 0 for z «0 was 
invoked. This allows substitution of 


into equation (2.22.0). 


When (2.22.1) is evaluated at pop =b it becomes 


H“(a,w) = yb[B(a,w) I (Yb) - Cla,w) K (yb)] 


or 
H (asw) - B(a,u) yb I (yb) 


Ον Yb K (Yb) 


(227) 
If we replace the coefficients in (2.15) and (2.16) by the ex- 
pressions given in (2.26) and (2.27) respectively, then it is apparent 


that 


η ω 
Hg, (pa) = <> yb K yb) κ] (yp) 


and 


+ 
Vp, (Prasu) = Blasu) I, (p) - Hk Qo) 


By evaluation of the preceding equations at p = b and simplifying, 


we get 
H (asw) K (yb) 
H* (b,o,u) * H; (b,o,u) = - em 
$1 os $4 23 yb K Y 


A (asw) > 
Yip, raul + Py αν ,0) edo rens 
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Subtraction of the latter from the former yields 


$ . + = n B B A yu 
Hy, (Paw) Y, (bras) dial Y, Bra) BR tb) 


It is seen that 


+ + Ν 
"oy (b 5 5 W) - Va, tP 55 W) 


iM(w) 
Vàm b(atk) 


by evaluating (2.12.5) at p = b and using the equality given by 


(2.19). Substitution of this result leads to 


iM) — H. b TNR {ρ 7 B(a,u) 
Yon b(atk) i - ,0,0) n ,0,0) yb K (yb) (2.28) 


Equation (2.28) is analogous to (2.25). These two equations are the 
result of matching the tangential electric intensity (-9 « z « o) and 
the magnetic intensity (z » 0) at the common boundaries between the 
regions into which the space was earlier divided. These equations are 
also combined to form one of the Wiener-Hopf equations by eliminating 


B(a,w) between them. The resulting Wiener-Hopf equation involving 


W'(a,u) is 
iy M(w - W {anu 
ze [kK (ya) - K (yb) ] = er + yb K (yb) 


x [Hp (baw) - Ya, (Peau) ] 7 va Kalk, (2,00) - 9s, (a,0,0)] 
(2.29) 
Since C(a,w) in this problem is unique, the two values--one 


given in expression (2.24) and the other given in (2.27) must be equi- 


valent. When these expressions are equated the results relate the 
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two variables oa) and Theses n 


; _ a Katya) I,(yb)  I,(va) 
W (asw) = BK ey (a,w)] -Bla, oR ng] - CER K ow) 


+ 
The terms W (a,w) and B(a,w) appearing in the above expression are 
eliminated by using expressions (2.28) and (2.29), giving a Wiener- 


Hopf equation involving the function He (asta) Š 


iy Μω) py o (Ya) - I(yb)] = Hato) 
verlo + k) yb K (vb) 


+ yb Lo (ve JH, (sos) - PCE 


- ya Lo Cra DH, C500) - V4, (asasu) (2.30) 


In the coaxial waveguide antenna we usually have b-a << 515 

to ensure TEM mode excitation at microwave frequencies. For this case 
the field intensities will be rapidly varying when compared to p „and 
the fields in the annular region will approach those of a parallel plate 
waveguide [17, p.329] of width (b-a). 

For the parallel plate waveguide it is well known from the pre- 
servation of symmetry that for a TEM mode excitation the current dis- 
tribution on the two parallel plates are equal in magnitude and oppo- 
Site in direction [21, p.126]. For excitation by En mode fields, 


the currents on the two parallel plates will be equal in magnitude and 


in the same direction for odd n, and in opposite directions for even n. 
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Up to this point no approximations have been made and equa- 


tions (2.29) and (2.30) are exact. It is recognized, though, that the 
form of these equations is such that they cannot be decomposed into 
two expressions, one analytic in the upper half plane and the other 
analytic in the lower half plane with an overlapping strip in which 
both functions are analytic. This is made apparent by observing that 
the coefficients multiplying the "minus" functions on the right hand 
Side of each equation are not identical. These same coefficients also 
contain a branch point and are therefore not of the form where an 
exact solution is possible [4, p.153]. 

We note, however, that if b-a << an then the field inten- 
Sities approach those of a parallel plate waveguide and if TEM mode 
excitation is impressed, the currents on the cylindrical waveguide 
conductors are approximately equal in magnitude and opposite in direc- 


tion. The mathematical expression of this condition applicable to the 


structure under consideration is 
2nb[H, (b,z,u) - H, (b,z,0)] = 2ma[H, (a,z,w0) - H, (a,z,w)] 
$1 $5 ve Po 


for z «0 . We now recall equations (2.11) and (2.7): 


ΠΠ 2.11 
Hg, (052.0) P * Hg, (Paz au) ( ) 
A _ Mw) okz (2.7) 
p 


Substitutions of these quantities in the above equations yields 





e 

>> = E 
>> -ω £&- 
<> > 





E 


PH, (baza) = 9, (bz,u) - "iol eT] 


= aLH (a,Z,w) (a,z,u) - Mu) „ikz 1 tor z 0 
3 


- Yo, 


Under the assumption that the approximation is very good, we 


can equate the two currents and obtain 


b[H, (b,Z,w) - v, (b,z,u)] = a[H, (a,z,u) - v, (a,z,u)] 
$1 $^ $4 $5 
for z<0 
When this equation is spatially Fourier transformed the results are 


precisely what is needed to allow decomposition of the two Wiener-Hopf 


equations (2.29) and (2.30), i.e., 


Ρ[Ηρ (Ῥναω] E Ve, (boa) ; ων = dg, (35050)] (2.31) 
when 


b-a << 22 (2.32) 


It is also noted that the expression (2.31) is analytic in the 
lower half plane given by τς Κο ΟΥ k Re(y 1 - (8/2) depending 


on the excitation. It is convenient to define S (a) such that 


S (a) = BER (basm) - ALLE = a[H, rae) a 
(2:93) 


To derive decomposable Wiener-Hopf equations we substitute 


(2.33) into (2.29) and (2.30). The resulting equations are: 
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iM(w) » μ΄ (α, ω) Ta 


> E (2.34) 
vématk) yal (va)[k (vb) - K (ya)] 
iM) Ha) a) 0 (2.35) 


Verlatk)  v^bK Gb)LI (vb) - 1, (2)] 


It is interesting to note that if the radius of curvature is 
large we may imagine the conductors in Figure 2-1 to be cut along the 
Z direction and flattened into parallel planes. Therefore it is 
reasonable to expect that equations (2.34) and (2.35) will become the 
Wiener-Hopf equations for two parallel planes in the limit as radii a 
and b approach infinity. The planes would be b-a distance apart and the 
axis of symmetry would be displaced a distance a=b >œ from the 
coordinate z axis. By taking the limit as a and b >œ and using the 
asymptotic foms for the modified Bessel function, equations (2.34) 


and (2.35), respectively, become 





iM(w) _ κου = m ES Ss 1509 
vén(a+ k) E προς ο πο) 

and 
iM(w) " μον [ - Yb O 
meer) = : sinh ra 


+ 
It can also be shown from the definition of Wlaw) and H (a,w) 
+ 
that for a parallel plate case H (aw) = - W (aw) . Therefore the 
equations (2.34) and (2.35) become identical in the limit as the radii 


approach infinity and each is the Wiener-Hopf equation for the 
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radiation from an open ended parallel plate waveguide with TEM mode 


excitation [4, p.107 or 21, p.128]. 

The solution of equations (2.34) and (2.35) for Ν΄ (α.ω) and 
H* (a,u) allows a complete characterization of the magnetic intensity 
in the spatial frequency domain. The foregoing work already contains 
the equations relating Ala,w) and D(a,w) to "MUR and αυ]. 
These coefficients are given by (2.26) and (2.23), respectively. Το 
obtain the remaining coefficients C(a,w) and D(a,w) in terms of 
ο ο) and H’ (aw) , we retum to equation (2.22.1) and evaluate 
it at p=a andb . When these results are considered with (2.20), 


(2.21) and the relation 


d = 
E Low (0,0,0) ] = 0 
dp 92 p=a,b 


we get 


tł 


W"(a,) = ya[B(a,w)I (va) - C(asw)K (ya)] 


H (aw) 


yb[Bla,w)I (yb) - Cla,w)K (yb) ] 


The solutions to these simultaneous equations yield the coefficients 
we seek, 


H^ (a,u)K (ya) W (o, u)K, (yb) 


] 
Blase) = [Sp = ove ava) = Trak) 


I 
0 0 0 


+ + (2.36) 
H (o,u)I (va) W (o,u)I (yb) 


a) = I aK a) Ta G5)! 


(2937) 
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To complete the list we add: 


+ 
Ala.) = 34 LO 2.26 
asu) = FR OT -- 


W” (o, w) 


D(a,w) = Τα] (yay (2.23) 


We return to the earlier discussion that certain values of 


ka,kb are excluded. Namely 


ka = Uon 
and 
kb = u 
om 
where Uns: is the nth ordered root of J(u) = 0 
; a _ 
and Um 15 the mth ordered root of J(u) -J (5 Y) zd 


There is an important mathematical reason for excluding these particu- 
lar values of ka and kb that shows up in the factorization procedure 
of Appendix D. The physical reason for their exclusion is that at 
these frequencies the magnetic intensity everywhere outside the struc- 
ture is zero. 

If ka- Uon ’ then the magnetic intensity in the region 
o <a ,2z< 0 must be zero. This value of ka corresponds to the 
cutoff wavelength for the E mode of a circular waveguide of radius 
a. Exactly at cutoff, the magnetic intensity for the Eon mode in 
the circular waveguide is zero [17, p.322]. 

Now considering that b-a «« 555 the structure has fields that 


behave approximately as those of a parallel plate waveguide. In this 





mos 
limit we have shown above that the field exterior to one of the plates 
is zero. By symmetry the field that is on the exterior side of the 
other plate must also vanish. By continuity it is concluded that the 
field is zero everywhere exterior to the waveguide channel. 
For the other case, kb = Um > we examine the asymptotic 


behavior of |u as moo 


jl 





Jo) - NC u) T cos(u - = - 


[ω]το 
ς 


or 


a T 
cos(— u - $) 
à b 4 


-JE (cos(u - = - cos ($ UNE 3) 


with a and b of the same order of magnitude as u >œ  . Through 


the trigonometric identity 


COS X - COS y = - 2 sin) sin = 


it is determined that the zeros are approximately located at 


2Tb 


2Tb 
lu | = {π-} (650) or m) , m-1,2,5,--- 


T 
> © or mo, 
for Uom 


The asymptotic values of u given above are recognized as 


om 
the asymptotic values of the cutoff frequencies as m +» for the 


: * 
E om modes of a circular waveguide of radius ar and the even 


ordered E om modes of a parallel plate waveguide, with a channel 
width of b-a , respectively. The second value of D also cor- 
responds to the asymptotic value of the even ordered E om modes of a 


+ . 
coaxial waveguide [22]. In the event b-a << a52, we may approximate 
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the structure as a circular waveguide of mean radius ae or a 
parallel plate waveguide, depending on the excitation. At the cutoff 
frequencies for such structures, the electromagnetic fields are iden- 
tically zero. Therefore we have excluded the use of such frequencies, 


Since they invoke no response. 


Solution of the Wiener-Hopf Equations 


The two Wiener-Hopf equations (2.34) and (2.35) are solvable in 
the conventional Wiener-Hopf fashion. The various terms are factored 
into the product of a "plus" and "minus" function, where required, and 
the resulting equation is decomposed into the sum of "plus" and "minus" 
functions. [The "plus" ("minus") functions are regular in the entire 
upper (lower) half plane.] Analytic continuation arguments are then 
invoked to yield the solution. 

To preserve the continuity of the work, Appendix D was devoted 


to the factorization [23] of the kernels 


x (a) x7 (a) 


(1, (ya)[K, (yb) - K (ya)]}7! = x(a) 


and 


{K (vb)EI (vb) - 1,2) F y(a) = y'(0) y (a) 


+ - + - 
In this section we will use the expressions x (a), x (a), y (a), y (a). 
Their specific functional form is found in Appendix D. 

We rewrite (2.34) and (2.35) using the newly defined symbols 


for the kernels that must be factored, 
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ole W (aw) * - (u) = 
Zark) 2, [x (a) x (a)} + S (a) = 0 (2.38) | 


7T (atk) A ly (a)y (a)l + S (a) - 0 (2539) 


Recall that the value of ης is given as 


When this value of e is substituted in the above equations, we can 
reduce them to 
+ + Ξ 
ἱΜ(ω)(α-Κ) N (asw) x (a) + > (a)(a-k) _ 0 
Vàn(o*k) x (a) ala + k) x (a) 


and 


iM(w)(o-k)  , Hita) y (a) , S" (a)(o-k) . g 


v2r(o+k) y (a) b(a* k) y (o) 


Examination of the region of analyticity of the terms in the 
above equations reveals that the middle term in each equation is 
analytic in the half plane t» -k, and the last term in each equa- 
tion is analytic in the half plane T < Κο . The first term in each 
equation has a pole at a=-k . If this pole were not present it 
would be analytic in the half plane tT < Κρ . By removal of the 


residue of this pole, we obtain the completely decomposed equations: 
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W (aw) x" (a) -į ᾱ. ΚΜίω) 


alatk) ~™ (atk) x (k) 
p [S (a) (a-k) m iM(w) [tk n 2k 1H 
x" (a) VOm(atk) x(a)  x'(k) 
H (a, w) y ME ) | Bo — 
b(atk) (otk) y 


2.40) 


= - (Sm (ok), iM) po-k , λα. 41) 


y (a) νσπίατκ) y (a) πμ 


In the above expression we have used the result from Appendix D, that 


X (o) ENSE) 
and 
y (a) » y (-o) 


The terms in (2.40) and (2.41) are arranged so that the left hand side 


is analytic in the half plane το 2 =k, and the right hand side is 

analytic in the half plane τρ « Κο or k perque) depend- 
ing on the type and frequency of the excitation. The important point 
in this discussion is that there is a strip of overlap in the regions 
of analyticity, that is, both sides of (2.40) and (2.41) are analytic 


in a common strip 


-k < T< Κρ or Κρ Re J1 - (8/k) 





AT. 

We therefore conclude that each side of the expressions (2.40) 
and (2.41) represents the analytic continuation of the other side. 
Since these expressions are regular throughout the complex plane they 
are an entire or integral function. 

If we can show that each side of (2.40) and (2.41) is bounded 
as |al + in the half plane for which that side of the expression 
is regular, then we can use Liouville S theorem to set each expression 
equal to a constant. Liouville's theorem states that the only bounded 
entire functions are constants. 

To determine the limiting behavior as |a| +œ for the terms 
in (2.40) and (2.41), we examine the edge conditions given in (2.14) 


as 


55 LOH, (0.24) ] a ze 


as z^ o at o= aorb . Taking a slightly more general case, the 


Fourier transforms 





Taz 
7" e ^^ dz 


and 


n laz 
Z a dz 


X 2 
83_—— O O——8 


A |- 


determine the limiting behavior. In these integrals n <0 to have 


absolute and uniform convergence of the improper integrals. 


E 
It is also well known that if |f(z)| «e ? as z^ then 


; : : + 
lim f(z) = lim iaF (a) = (lim ia)( lim F (a)) 
Ζ 0t a> © (1 > 00 a> © 
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KZ 
τ τ» -K, where a=otit . Similarly, if |f(z)| « e e as 


z > -œ then 
lim f(z) = (lim ia) (lim F (a)) 
z >07 Q -> =% Q > -%0 
for t<k, . Thus the limiting behavior as |a| >% in the region 
of regularity of the spatial transforms is determined by the behavior 
of the space domain function in the limit as the origin is approached. 


Using the gamma function [18, p.183] 


oo 


"ms | Er σα 


and contour integration, we find 


ae = οἶαΖ posee .] r(n*l) A z(n+1) 
E Jo" 
and 
0 i z(n*1) 
ES | Meia c 1 T(n+1) ec 
ver τ αν 


OO 


with -1 <n <0 in both integrals. The lower extremity on the 
range of n is set at -1 to ensure the existence of the gamma 
function. 


From the edge condition we have n-- 7 and determine 


1/2 


Was) ν a 


and 


Hsu) v «^ as jal >= with 1 > -k, - 





> 4 p 
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From the definition of S (a), it can be seen that it is the 


spatial transform on negative values of z of the current normal to 
the edge at z = 0 on either the inner or outer conductor, exclusive 
of that current contributed by the incident field. The magnetic in- 
tensity is continuous at all values of p for z=0, thus the 
current in each conductor must approach zero as z>0 . This is 
apparent from the fact that the current in each conductor is propor- 
tional to the discontinuity in the magnetic intensity across the con- 
ductor. The discontinuity disappears at z= 0 . Because the total 
current approaches zero as z>0 , the quantity which transforms to 
S (a) must approach that of a TEM mode field intensity traveling in 
the minus z direction. This quantity should just cancel the incident 


mode field intensity at z=0 . Thus 


¿mal \a.2,0) E Up z NETRE D 9g, Vo z,u)] 


=> M(u) eus as z>0 


The transform of this quantity is 


0 
: ο. -ikz az, _ -iM(w) y 1 
25S (a) 'ν p [πω ο a dz = cU Lh 


TT 


=O 


From this result we conclude that S (a) a (E as la] >= with 


m< tko ; 
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Collecting some results from Appendix D, we write 


x (a) A [LE 
and 
y(9 * la] V^ as Jal >e with T> =k, 
Also, 
x (a) v PE 
and 
1/2 


y la) v lal as |a| » with t< k, 


An examination of (2.40) and (2.41) in the light of the 
asymptotic behavior of the terms shows that each side of these ex- 
pressions goes to zero as |a| >œ in their half plane of regularity. 
Liouville's theorem is therefore applicable and each side of (2.40) 


and (2.41) must be identically zero. The quantities of interest are 


Wlaw) χα) a E kM) -o 
alak) ο i JZ (ak) X'(k) - 
and 
H (aw *a) . + [2 _k Mu) -0 
τος ee 
Therefore 
a iJ ee (2.44) 
x (a) x (k) 
H^ (au) = j 2 kb Mw) (2.45) 


7 y'(o) y (k) 





sr 
and the coefficients A(a,w), Bla,w), C(a,0) and D(a,u) have been 
determined, see equations (2.23), (2.26), (2.36), and (2.37). The 


spatial domain field intensities are now written: 


ο .P k M(w) Ky (yo) 
Η, (ριαιω) = - i |Ξ -------- ------------ (2.46) 
$1 ? YK (vb) y (k) y (a) 


= 1, (yp) 
E ums 2 um EE Fc. 2.47 
esa μ' vI (va) x (k) x lo) P" 


and 
. [2 KM(w) | | 
Vy, (o. * i JE PL ο πο. 
bo x Y I, Y Ko ya LA ya Ka Y 


I(yp)K,(va) + Kılyp)I (va) | Ij (vo)K (vb) * K (vo) Lo (yb) 


} 
y (k) y (a) x (k) x (a) 


(2.48) 
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3. Frequency Domain Description of Fields Interior to the Coaxial 
Waveguide 


In this chapter we shall derive the frequency domain fields 
interior to the coaxial waveguide. We will completely characterize 


the TEM mode as well as the EN (q=1,2,3,---) mode fields. As pre- 


q 
viously stated no other modes are allowed due to the symmetry of the 
structure. 


Let us begin by recalling equation (2.11) 
H 9 9 = η + 9 9 
p P Z ,u) 4 Ve, Vo Ζω) 


The incident field is known to be 


yl _ Mw) eikz 
$ p 
Our task is to find Y, (p,z,u) . From the definition of the spatial 
2 
transform we write 


Py (PZW) = ae | Y, (p,a,u) e” ἀα (3.1) 
2 ven ὁ 2 
Before attempting this integration we will examine the integrand 
and all singularities will be located. To do this we make use of the 
relations [20, p.375] between the modified and ordinary Bessel func- 


tions that follow: 


- ; vmi 5 mi - 

1 (2) =e J (z e ); (-n <argz < 3) 
| vri Ti 

K, (2) = 7 mi e? Tr e? ); (-"T « arg z ELS 
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We also make use of the analytic continuations [24, p.80] 
mri. . muri 
e e NEA 


mmi, | .-mvrri . Sin mvm 
K (z e )=e Κ (2) Ver I (z) 


By these formulas it can be shown that 


ΠΠ Ta) uti» 7 — πρίν ΒΥ. (γα) ZU man ya) J 


where Υ = -i Jk*- E = -iy'. 


It is clear that expression (2.48) does not have a branch point 
at a= +k , Since we may replace y by -y and not change the value 
of the expression. However, (2.48) has a branch singularity at 
a = -k due to the presence of x" (a) and y (a) . From the above, 
we conclude that the only singularities of ne which exist in 


the upper half plane are poles, located at a=k and a= a 


q 
(q=1,2,3,---) with 
EZ . 12 EN? 
Og = /k - Um 1 Yoo k 
where Y oq is the qth ordered root of 
Jo(v'b)YoG'a) - Jo (y'a) Y (y'b) - 0 e 


The oe determined by (3.2) are the eigenvalues determining the cut 
off frequencies of the E og modes of a coaxial waveguide [22]. 
To evaluate the integral (3.1), we integrate along the closed 


contour in the upper half plane as shown in Figure 3-1. On this 
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Q plane 


» 
> 
>! 
n 
? 
EU 


Figure 3-1. Integration Contour 
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contour the integral converges for z « 0 and the contribution along 


the semicircular path at |a| +» is easily seen to be zero by re- 
placing the modified Bessel functions in 9g, oro) by their 
asymptotic values and using Jordan's lenma to evaluate the resultant 
integral over the portion of the contour at |α| » œ. 


By the residue theorem, we have 


Yo (p,Z,w) = 2111 ) residues of (— v, (o ,a,u)e ! 97] 
2 ym ϱρ 
at the poles interior to contour 


of Figure 3-1 ο] 
The residue of the pole at k is 


Hin(o-k) iE Ge y 


|; I_(yb)K, (ya) - ΓΚ τ᾽ 
I,(ye)K (ya) +k, (ye )I (va) | I, (ye)K (yb) + Ky (yo) I (yb) 
y (k) y (a) x (k) x (o) 
= ik Muje 7 MEE) 
a> k Yatk 
[y OIA an yar) - OT ΧΡ 1π γῦ 4 το) 
mun ο... ο ο ο AA YP 


gn b/a 


x 


-ikz 
er? - RT? E (3.4) 


where limiting forms of the Bessel functions have been used to charac- 


terize their behavior for small arguments. 
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The residue of the pole at ας is 


å -iaz 
lim (ara, ) ikMu) e 


a> my{~ ald (y'b)Y (y'a) - J (y'a)y (y'b))]} 


IH (ve)K (va) * Ky(ve)Jglva) — T4(ve)K, (vb) * Ky(vo) Jo (yb) 
y () y (a) X (k) x (a) 


Using the L'Höpital rule this becomes 


-iz 
lim [- kM(w)e ] 
α” αι B 


2 + k) + N E k) UT 
E Ey! (- 7) (J (v 'b) Y (y'a) Ν J, (y'a) Y, (y'b) )] 


(3.5) 
Observe that y'= kê- a and Yoo = yké- a. . By invoking relation 
(3.2) the derivative in the denominator of the above relation is readily 
reduced to 
d 


E bre 5JL9 (v ' b) Y, (v a) - AN 


Γ Ζ υ'[{π{0ρ(ν'Ὁ)Υρ(γ'α) AA 'b2)] 


t 


N ΣΥ [ττ][ο(ζη(γ'ὈΥρίγ'α) ~ doly'a)¥oty'b)) 


+ 


a(J (y'b)Y4(y'a) - J4(y'a) vo (y'b))] 
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Putting this value in (3.5) yields the residue 


| -izJk^- (ντο 
τα ολλ. Ya rg] 


oq 
E aay 


J (y! b) 
0 ‘og 
1 Ye Yoa? 


UE RIA AR o) 


A ad Yo (Yoq a) - O A AAA) 


x 


(3.6) 
A, ας b) 
The relation ο - ο was also used to simplify (3.6). This 


o'Yoq 
relationship follows from .. 


ve (0,2,0) for z «0 is found by substituting (3.4) and (3.6) 
2 


into (3.3). The result is 


Ve, Vo 2.) -M(u) Ree - ) MR f, (o) Yog 


R - gpz [D OT? - OT (3.7) 
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M αΚ 


| Wk 2 A oq 


| Jo Y esp) 
| y UKY SKE- rE) Toog?) xk) xE A) 


ee reg) 


(3.8) 


fate) = δη (ραΡ)Υ ο( 1ος) - Υ1{οαβ) Yoga? (3.9) 


The magnetic intensity for the interior of the feedline is seen to be 


Μ(ω) a KZ 


Ha, o 2,0) - 5 


+ ὑρρίρ»Ζ»ω) αμ ς ο (3. 10) 
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4. The Radiated Electromagnetic Field 


To find the magnetic intensity in region 1, we must take the 
inverse spatial Fourier transform of the spatial frequency domain mag- 


netic intensity found in Chapter 2. The integral to be evaluated is 


k M(u) Ky (yp)e 1% 


o (4.1) 
yK (Yb) y'(k) y (a) 


P 
Hy (0.2.4) z T | 


for p >b . This equation was obtained from equation (2.46). 

Examination of the integrand of (4.1) using the analytic con- 
tinuations of the Bessel functions and the value of y on both sides 
of the branch cut reveals that the integrand has branch points at 
a=tk. Any attempt to evaluate the above integral in closed form by 
contour integration would have to account for integration along a path 
alongside the branch cut. If the contour is closed in the upper half 
plane (for z< 0) then the contribution from the path alongside the 
branch cut in the upper half plane must be calculated. If the contour 
is closed in the lower half plane (z > 0) , then the contribution from 
the path alongside the branch cut in the lower half plane must be 
calculated. In the special case of z = 0 the integration around the 
branch cut in the lower half plane resulting from closing the contour 
in a negative sense equals the integration around the branch cut in 
the upper half plane which results from closing the contour in a posi- 
tive sense. There are no poles located inside the contours under 


consideration. 
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The integrations around the branch cuts are quite complicated 


and are not evaluated in closed form. Since our primary interest is 
in the radiated fields (i.e., those with 1/r behavior),it is sufficient 


for our purposes to approximate the value of the integral (4.1) for 


r>0o°5 


In the limit as r ^» , we replace the Bessel function 
dependent on p by itS asymptotic value for large arguments [20, 


p.378]. Using Figure 4-1, we readily have 


p= r cos x and z=rsin x 


where |x| < 1/2 . Using the asymptotic value 


manm cp 


the integral (4.1) becomes 


H, (r,0,w) = da 


-iM(o)k | e" Yr cos x * iar sin x) 
$1 


2s y (k) YK (vyb)y (o) LArr cos x] 
(4.2) 


as ro, p»b 


The variable @ has been introduced at this point since this 
is the variable in which we wish to have the final expressions. It 


is obvious from Figure 4-] that 





BG ie 


| (7,0, 4) 


(p.d, O) 


Figure 4-1. Coordinate Axes 
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Define a change of variable by introducing the transformation 


a= k sin z (4.3) 


z is the complex variable xtiy . Under this change, the function Y 
becomes 

y = J/k*sin?z - k^ - -ik cos z (4.3.1) 

It is necessary to specify precisely the value of z , since it 


1s given by the multivalued function 
Z = arc sin(y) 


which has infinitely many branches. The precise specification will 
allow a one to one mapping from the complex a plane to the complex 


z plane. Since 


E zu e 17 
ru sin Z = 21 

we may write the equation 
o1? 2] 8αὶΖ.9 - 0 


This quadratic is readily solved for z 


z = -i [i c [1- (Θ΄ 


In this solution, the principal branch of the logarithm function has 
been used and the plus sign has been used with the square root function 


[25]. By the branch previously defined for y , the above equation 
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simplifies to 


z = -i &n[i BI (4.4) 


The equation just developed is useful in determining how the two- 
sheeted a plane maps into the z plane. 

Under the transformation (4.4) the numbered regions of Figure 
4-2 map into the correspondingly numbered regions of Figure 4-3. The 
path of integration for the integral of (4.2), depicted in Figure 4-2, 
maps into the path shown in Figure 4-3. 

As examples consider two points on the path of integration in 
the a plane 

1 o» -9 with [τ] » 0) 


2) o >+% with |τ| » 05 
For case 1: 


2 
= Tim ((-Jo] + i]t] )*- κ᾿] - lim Clo](1 - +5) - ilel] 
g >% g>% 2|o| 
| «| »0* | t|>0" 


and 


: + 
Z= -1 DAR anpitelol +41 113. Γ zs RS q =-7 + jo 


[50 
For case 2: 


2 
y= Tim [{|σ|-1|τ[]ξ- #232 = 1im Elo (1 - Sy) -iltl] 
l|» |o |> 2|o] 
[τ] 0” [τ]-- 07 


and 





300. 





Re y >O 


ES Im y >0 
[ ] Iny «0 


Yj). 






Uy 


Rey <O 


Figure 4-2. Integration path in a plane 
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z plane 


Rumes απ op Türe ----- E " 


"| s bp r 


Figure 4-3. Integration path in z plane 





Be: 

2 = -i Tim angi Le Haltn] = 4 timo anti lods 
|σ |< | o | > 
αυ. 


Upon substituting (4.3) and (4.3.1) into (4.2), the integral 


under examination becomes 


M(w) k 


H NOS = 
ον F VE y (k) 


< | - exp(ikr(cos x cos z - sin x sin z)) m (4.5) 
2 y (k sin z)K A -ikb cos z) V-irk cos x cos z 


The integral is in the form amenable to the saddle point method of 
integration, since r >% 


Define the z dependent quantity in the exponent as 
g(z) = i(cos x cos z - sin x sin z) 
- i cos(z * x) 
It follows that 
d er 
ο οσο τρις x) 


The saddle points z, are the points at which g(z) is a maximum. 


Therefore from the above we know that 


Zo = NT =- X ; Ποῦ, ο. ο 





ος 


Concentrating on the point z_ = -x , we note that g(z,) - d 


S 
Since a Steepest descent path is a path along which the imaginary part 
of g(z) is constant, the steepest descent path through ο ει is 
given by 


Im g(z) = 1 


Since g(z) = i cos(xtiy + x) = i(cos(x+x)ch y - i sin(x+x)sh y) , 


the steepest descent path through EE is given by 


cos(xty)ch y = 1 


This path is shown in Figure 4-4. We are assured that this is a path 
of steepest descent rather than one of steepest ascent, since the 
Re g(z) monotonically decreases along this path as we move away from 
the saddle point. 

By the residue theorem, we deform the path of integration of 
Figure 4-3 into the steepest descent path. The major contributions to 
an integral along a path of steepest descent come from points in the 
neighborhood of the saddle point. The slowly varying terms of the 
integral are replaced by their value at the saddle point and the func- 
tion g(z) is replaced by the leading terms of its Taylor series about 


the saddle point. 


yz) = i- Hz) 


ikr akr 2 
M(w)k e' | exp(-1 SAZtx)") dz 
H¿ (r,0,0) = 3n T 
| Jan y (k)y (-k sin y)v-irk cos x K (-ikb cos x) 
(4.5) 
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ly 


Steepest descent path 


Figure 4-4. 
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Along the steepest descent path 


e 1 
CM cos (x+y) 


dy _ ,„ _tan(x+y) 
dx es y 


and 


As we approach the saddle point in the manner x>-x and 


y Zu , the slope is 


As we approach the saddle point in the manner x -x and 


y>0 , the slope is 


lim Y = -] 
X -X 
y >0- 


Therefore the angle that the SDP makes with the real axis at the 


saddlepoint is either 


3m T 
rum. 


Examination of Figure 4-4 shows that at the saddle point the SDP and 
the negative real axis form the angle 37/4 in the second quadrant 
and -7/4 in the third quadrant. This observation permits us to 


express ' 
| | u| ΗΡΙ ; ue SDP and y >O | 
Zot == 


| n σα, π/8 


> ue SDP and y-0 | 
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It is also observed that 
u- eu e 7/4 at the saddle point 
Since equal contributions to the integral of (4.5) are received as we 
move away from the saddle point in the two directions of increasing u, 
we write the integral which appears in (4.5) as 
ο kr ΝΠ 

Sul? -iF 

e 


| exp(-i S (z6)*)dz ὦ e d| ul (4.6) 
SDP 


OO 


Using the well known integral 


00 


e*9" dx= T 


fy 


the value of (4.6) is seen to be 


kr, 2 SET IR 
- Jul - i= E 
2 f(e ἘΠῚ Educ - ΠῚ (4.7) 


0 


Equation (4.5) becomes 


πα 1 
Hy (15850) = Mw) S— {— (4.8) 
1 y (k)y (-k sin x)K (-ikb COS x)cos x 


To simplify the above expression we use the identity relating modified 


Bessel functions with ordinary Bessel functions: 


IRR 
K C7 ikb cos y) = i > Ho (kb cos x) 





2 s 
From Appendix D, the expression 


] 
-} ΠΗ _ 
y(a) | > HS (Kb) [J (Kb) κ] 
is formed, where K = tiy . In expression (4.8) we have 


i z ut (ko cos x) = i HD (κ sin 0) 


If in the expression yla) , K is set equal to k sin 0 and a is 
2 


-K2 = 


set equa] to k cos 0 , which follows from the fact that y 


αξ - k* , then 


y(k cos 8) = y (k cos 0) y (k cos 9) 
-] 
| T HET) (kp sin o) [J (kb sin 0) - J (ka sin 1] 


When this equation is solved for i Z aC (o sin 8) and the result 


substituted into (4.8), the far zone magnetic intensity is seen to be 


alk J (kb sin 0) - J_(ka sin 0) wae s 0) 
om ws I TA | 


(4.10) 





The relation (4.10) has the restriction that p»b. In the 
region z <0 the equation (4.10) is only valid in the region p» b; 


however, for z » 0 , the H can be analytically continued from 


$1 
region ] to 2 and from region 2 to 3. Therefore the restriction re- 
quiring p > b can be removed in the region z » 0 and equation 
(4.10) becomes an expression for the magnetic intensity for all values 


of p in the far zone (r > ~) exterior to the feedline. 
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It is recognized that equation (4.10), except for the factor 
y (k cos 9) 


+ 
y (k) 
the incident field is the total aperture field and that no contribution 


, exactly equals the results obtained from assuming that 


is made to the radiated fields by currents which flow on the waveguide 
structure. 
+ 

The new field factor contribution al SEE to the coaxial 
waveguide antenna first obtained in this we the field pat- 
tem obtained by the conventional solution to the problem of radiation 
from this type of antenna [ 19,p.302]. It will be seen later that this 
factor has a significant effect at high frequencies and modifies the 
field pattem to direct the radiation in a more forward direction. 

By equation (2.11) we expressed the magnetic intensity in 
region 2 for all z as the sum of the incident field and another term 
which represented the reflected and higher order mode field intensities. 
By the radiation condition, we know that the incident field cannot 
exist in region 2 for z>0. 

To show that the incident field intensity is canceled in region 


2 for z » O0 , we evaluate the integral 


eco 


Y, (psZ,w) = —| va (0,0,0) e 
Po vem 4 ορ 


gH Uu. 250 (4.11) 
by using the contour of Figure 4-5, va (p,a,w) is given by (2.48). 

2 
The value of the contour integral is calculated by the residue theorem. 
The value of the integral (4.11) is the negative of the contribution from 


the path around the branch cut and -2mi times the summation of the 





=) oe 


+k 
TH en 
RV oy 





Figure 4-5. Integration contour 


«— C 





S94. 
residues at the poles of the integrand located on the interior of the 


contour c. The contribution along the semicircular arc at Jal >œ is 
zero. 

The integral (4.11) is not evaluated in closed form due to the 
complex integration on the path alongside the branch cut. However, we 
shall show that the contribution to the value of the integral from the 
circular path around the branch point (k) is just enough to cancel the 
incident field intensity. 


16 


3 


On the circular path around the branch point œa = -k + ee 


e>0 and 0<6<2n. The integral of interest becomes 


tim |- P e 192 2 
€ 0 von 
= 


lim Zu YO: 


: 27 1 1. d 

Í -ikM(u) eif | psc “πιῶ 
; F F + 

ε- 0] τ an = 0 y (k) y (a) x (k) x (a) us 


; 2T | 

= Jim d zik MG) eK? eT E 1 
AL a 

e>0 7 £n c y a 0 yo | 


. ο Μ(ω) 
D p an È 


> an ET o 


We have made use of the limiting forms of Bessel functions for small 


arguments and the following: 





EJ. 


x (-k) = x(k) 5 χ᾽ χα) 


x(k) 


t 
I 
er 

i 


N 


y(k) 5 y'(k) y (K) - y(k) 


To see that this is just enough to cancel the incident magnetic inten- 


Sity, note that 


which is easily found from the definition of y(k) and x(k) given 
in Chapter 2 and the limiting forms of Bessel functions for smal] 


arguments. 
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5. Time Behavior of the Radiated and Reflected Electromagnetic Fields 


The time behavior of the field intensities 1S found by comput- 
ing the inverse time transform of the frequency domain expressions. To 
find the time behavior of the radiated magnetic intensity emanating 
from the waveguide antenna, the inverse transform of equation (4.10) 
must be obtained. Because of the complicated functional behavior of 


+ 
the factor Uk cos 9) the computation of the integral 





y (k) 
ikr J (kb sin 6)-J (ka sin 0) 
E a] 2 O O 
ο τι (2 gs —— — 
+ : 
y k 3 8) e iot au (5.1) 
y (k 


is extremely difficult. 


For harmonic time dependence the exciting voltage V(t) is 


iw t 
Re(Ae ° ) (see equation (2.7)). Then 


M(w) = R zm fe A δ(ω-ω.)] 
(w) = Re[ VZn u, In b/a ω- ω, 


where S(x) is the Dirac delta function. Therefore, the case of 


harmonic time dependence (5.1) reduces to 





nn [en N 
Ho en b/a y (2 r sin 6 


Wb l wa 
x I nn 0) - ο πε π 1| (5.2) 





LI. 
The value of k is w/c as given by (2.4.3). 


Tow 
The factors of equation (5.2) except y tg. cos 8) are well 


T 
y (C2) 
understood, and are not in need of further simplification. 


From Appendix D we transpose the function 


K 9 
| + ———- . (at 
y (k cos 0) 5 ! [et] | (ΑΡ) Κ(1. cos 6) 180k sin 8 
- k cos 6 e 
y (k) m=1 | 1+ 05 9 
i | oml 
att | Am : 
Vee ie k(=5=) (cos 6- 1)(1-C+i 5) 
x (exp f KG) (u) an) ke NM. e 127 2 
EE 14K cos 8 
Vk? - ul 


atb ) atb 2b 
ep (E cos OL -2n(k (5) cos 0) +35 2n(2 cos 0) 


X 








à ib) HM ma) (5.3) 


In this formula we have used the fact that if a=k cos 0 , then 


8 -ik sin 8 + k(cos 6 - 1)4 _ 


+ (a-k) Πα. .; 
nern = Alk eins kteos = 1, * la] > -i0 


e 


using the principal branch of the logarithm function. Additionally, 


(ΖΗ: 2 | 
K — — ] m — περι. 
di Zi ΠΕ emp 
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It is also known from Appendix D that y (a) " [α| 46 


as 
a9. From this asymptotic behavior we see for very high frequen- 


cies (à << b-a) that 


+ 
im ο eos ie T2 (5.4) 


+ 
k> œ y (k) 


In Chapter 2, there was imposed a restriction that kb,ka << 1 
for non-harmonic time dependence. In keeping with this restriction, 
equation (5.4) should only be used for the case of harmonic time 


dependent excitation. For values of ka,kb << I the value of 


y,(k cos @) 

EU) y is more complicated but instructive to determine. To 
+ 

determine the value of this field factor in the low frequency limit 


(à >> b-a) we use the series expansion of the infinite product 


en - cos 0) 


"περ 


= ] - (33) (1- COS 0) 5 IM L T olka)? 


m=] 


The infinite series is the difference of two divergent series but is 
itself convergent [26, p.33]. This is apparent since 


] E BE cb... 1l. e 
atb m 


o 1 
[anl Bo (m -q Anm 
as mo and the infinite series } L is known to be absolutely 
m=] m 
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convergent. The following integral 


k 
(2) V kt uc 
K`’ (u) £n o du 
D een 


ὁ 
Vkĉ- " 


] + 


also appears in en. This integral will now be shown to be 
convergent and its OO value will be determined in the limit 
kb,ka < ]. 

Upon substitution of the value of Kl?) (u) and simplifying, 


the above integral becomes 


f πμ ος 
$30" ΠΠ cos 0 +y¥1- (if 


Introducing a change of variable, ub-z , we get 


Do dete 1 J| onfi + l - cos 6 la 
sot 32(z) + γ5(2) cos 0+ f= (ie 
(5.5) 
To examine the behavior of (5.5) for small values of z , we 
consider 


€ 
] 2 ] l - cos 0 
A 


coso+y1+ (£^ 








€ 
E. ] 1 
This integral converges as | iz mom dz converges. 
---- {η z 
0 2 
T 
Now 
€ 2€ 
A] dz npa 4 
I+z an Ζ ὁ 2 {η (7) 
TT 
€ 
= ne d[ - 1 ] =. τό ] - 
4 Qn z 4 nz 
0 0 
This quantity is seen to be convergent for all ¢€ # 1. We can also 


remove the principal value sign and change the lower limit to zero 
without affecting the value of the integration. 


At the other extreme for large values of z , we consider 


1 f 2 1 l - cos 6 
— Ir — Da 0 a 
a TZ TD Gr o 5 


2 
cos 8- if (25) = 


for M >> 1 . The asymptotic value 
i(z-4) : 
ADO). (e [ιν ο] for z» 


is now substituted with the result 


- [1 - —— ] anfi mso IE 
[1+ 0(2)] E 
cos 8- i (cg) - ] 
en 1 Ἰ - cos 8 
"| 0(7) enf: νι: 
1 cos 0 - 1 (2 - ] 
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The logarithmic function £n(l*x) can be expanded 


2 
£n(l+x) = x - =~ + στ - | rer 


When this expansion is considered along with the fact that 


1 + Ἰ - cos 6 A 
Z 


σης 
cos 9 - 1 ui) - 1 
as z>% (with z real and |A| finite), then 
ai = 14; - 
cos 6 η. 
1 1, ‚A 
LY 0(7) (7) dz 
M 
This integral converges. 
The integrand of (5.5) has an integrable singularity at z=kb 
when 6= 7/2 . By the above we conclude that, for all values of 9 
and k the integral of (5.5) is convergent. 
y,(k cos 0) 
Since our interest is in the value of πε for 
+ 
ka,kb << 1 , we next endeavor to find the value of (5.5) in the limit 
as Κυ” 0. We will use this limiting value as an approximation for 


the integral appearing in equation (5.3). 





n - 


00 


lim 1| m- £ (71 40 (19 52 ) 62 
kb >0 Jo z) + Y.(2) A A V] _ (27) 
kb 
= lim :| [1 - > 5 £n (+ —— 3 dz 
kb +0 392 * o C2) cos e + / 1 - (E 
+ Tim τ [1 - EAM MÁS ο _ 2 Je 
kb >0 b o2) tv oz) Cuts a I ES cla 
kb 
= 0 
Note that the quantity 
lim [gn(1 + la TEE )] = 0 


kb >0 2 


COS 0 Jn - 1 
From the above work we conclude that as kb goes to zero, the integral 


(5.5) goes to zero. Thus, for low frequencies, we approximate the 


function k 
exp K2) (u) gn N UL. 
δ + n k cos 60 


Jk- u 


du 


appearing in equation (5.3) by unity. 
If in (5.3) we make use of the limit that x in x x as x90, 


+ 
then the field factor κος ο. can be approximated as 
y (k) 
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+ oO 
vases) ~ (1 = (5350 - COS o] D a ae 0(ka)?) 
= atb mm 
y (k) mI dol + 
x exp | i L sin a + Γκ [οί] - cos 0)]- £b cos 9 àn(2 cos 0)- £n 23) 
for ka,kb << 1 (5.6) 


The time behavior of the far zone radiated magnetic intensity 
for a non-harmonic time dependent incident field intensity with 
ka,kb << 1 is given by (5.1) when equation (5.6) is used for 
Y (k cos 0). In the far zone the value of r in equation (5.1) is 
Es E. larger than b or an To the approximation used, the 
exponential function in (5.6) is replaced by one and for that reason 
it does not contribute to the phase factors in integral (5.1). In 
addition the term of (5.6) containing the infinite sum can be dropped. 
This is apparent since ka,kb << 1 and the Fourier transform of this term 


will be of order 1/¢, which will make its contribution very much 


smaller than the contribution from the remaining term. In view of this 


[J (kb sin 6) - J (ka sin 0)] do 


for ka,kb «« ] (5.7) 
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When the incident magnetic intensity is harmonically time 


dependent, then the value of (5.7) is easily obtained. In this case 





M(w) = em [o — δίω-ω.) 
{η b/a 


and for low frequencies (X >> b-a) the radiated fields are 


Em cos (w (t- =) 


ο ΠΝ 
H¿(r,0,t) ~ A uU. ο ος να, J (—— sin 0) - Jose sın 6) ] 


0 


Ho 
E,(r,0,t) = [Fen fron (5.8) 


For non-harmonic time dependence the radiated field can be 
found by direct integration of (5.7) or by a time domain convolution. 
The integrand of equation (5.7) can be considered as the product of 
the Fourier transforms of two time domain functions. Therefore this 
integral can be evaluated by the convolution theorem and the impulse 
response is readily identified from the result. 

In applying the convolution theorem, use is made of the fol- 


lowing well known integral [27]: 


ve/n 
FR OE sz 
a vz -t 
B. | I.lazje te da = co : |t] = 
Br 
0 . . otherwise (5.8.1) 


x 


Upon evaluation of equation (5.7) by the convolution theorem we get 
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the result which follows: 


t 
H,(r,0,t) | M(t'-x) f(x) dt (5.9) 
where f(t) is 
0 s a 2 sin 0 
ΒΞ ‚2 sino<ı<- sino 
JR sin 8)" - τ 
ss ] ] A dcm 
f(t) = pam c 1 - ee Ok ED 
(b si ay? ED m 2 2 E 
sin - T (on 0) - t 
Ben ο... : = sin DU <2sin 0 
VE sin as T 
C 
| 0 ; T>—Sin 0 
(5.9.1) 
and 
ee ee 
t =t z 


The response to an impulse of voltage excitation is therefore 


Ea 
Hg (O9) 7 Ju ae sin 6 An b/a 50) 


From the character of the impulse response we see that at an 
observation point in the far zone, the radiation lasts a total time of 


(2b/c)sin 8. By ignoring for the moment the frequency restriction of 
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ka,kb << 1 and allowing an excitation with zero rise time, a step 
function of voltage , then by integration of the impulse response we 


obtain the step response, 


eo ] 
H¿(r,0,t) = wen’ ^ 
o mr sin 86 £n = 
= b de 
TT . =| τε b . : 
+ ποπ 5 =- — c c 
7 sın Ben ; E SINSO < T < sin 6 
rn A A ορ 
b sin 9 asino ? C En — C un 
i . =| τς a. b. 
. AT I = = 
7° SIN p-cin > ea” κ τς c 5n O 
τ | 
0 ; T275sin® | 
(5.11) 
r 
where Ben 


It is advantageous to have the step response for numerical cal- 
culations since its magnitude is finite at all points and the far zone 
fields can be found by the convolution of the time derivative of the 
excitation and the step response. This procedure does not violate the 
low frequency restriction. 

Because of the low frequency restriction the input excitation 
should not contain any high frequency components. To render the high 
frequency components in the excitation negligible, the rise time of 


the excitation pulses must be long compared with the time b/c . If 
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the input is of the form [1 - eR] u(t) then the time constant T 
must be chosen so that T >> b/c . This will preclude any violation of 
the low frequency restriction. 


For such inputs, equation (5.9) contains convolution integrals 


of the type 
C e 
l-e ll b D 
— dr, with - Su oirn z SİN 0 


Dr ee mr 2 
= ie be sin 8) - « 


and t< t 


t’ - |t-1| t' „ei 





| (l= es! Dan. | are) 
© sin oC sin 6) -? sine (2 sin 9)*- τ 
E | i dI 
t / 
mesi eSI) | e Sn BETTE 


- 2 sin 0 JE sin o)" -« 
C C 


0«|t'| «P sine 


Since T »» b/c we neglect the second term and conclude that 
the response to an input of finite rise time with time constant 


T >> b/c approximately equals the step response multiplied by 
i 
(1 E =: fly 
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Excitation 
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Figure 5-1. Radiated magnetic intensity of the coaxial waveguide 
antenna for (a) step function and (b) exponential excitation. 
Ordinate: Hy - Ht olg 
Time constant for (b) is cT/(b sin 6) = .1 
Antenna Dimensions: a .8b 


r £n > sin 6) (volts) 
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A basic assumption used in Chapter 2 is that the total current 


in the inner conductor is equal and opposite to that in the outer con- 
ductor. This assumption implies that the contribution to the radiated 
fields from all parts of the structure outside of the aperture are neg- 
ligible. Examination of the impulse response shows that the radiation 
appears to originate from the edges of the open ended structure 
located in the aperture plane. The outer edge gives a positive con- 
tribution and the inner edge gives a negative contribution. The 
radiation starts when the incident current reaches the edge and lasts 
a time proportional to the time it takes light to traverse the diameter 
of the conductor. The difference in sign of the contributions is due 
to the difference in sign of the current on the conductors. 

We now turn to the problem of determining the time behavior of 
the field intensities for a<p<b; z<0Q. The magnetic intensity 


in this region is given by 


i ikz -1kz 
uz t) - f eztot | 
$^ Vr p p 
£ N (Q1 
ιτ, 4 [do (5.12) 
q=] q q 
-jw t 


For harmonic time dependent excitation voltage Re(Ae : Va 


M(w) = eps Je x 8 (w-w 1) 
a 





π. 


and (5.12) becomes: 





Ον ΝΟ (5.13) 


where βίως), R (wo) and P are the functions (3.7), (3.8) and 
(3.9) evaluated at w= ως » İf applicable. ως Should not corres- 
pond to the cut off frequency of one of the Een modes of the coaxial 
waveguide, since our formulation is not valid for these frequencies. 
We also note that as De R(w,) + 0 , thus the TEM mode reflected 
field tends to vanish for very high frequency excitations. If 

ka,kb << 1 , it is obvious that for large values of z the higher 
order Eon modes are negliglble since they are evanescent in this 
frequency range. Since the time behavior of the incident field is 


given, we shall only be concerned with evaluating the integral for 


the reflected TEM mode fields or 


00 


1 Mw) ,.-ikz- iwt 
H το ο) προ | —— Re dw 


+ -2 + = 
SAH ="? (5.14) 
Qn a 


From the values of y (k) and x (k) determined in Appendix D, it 


can be adduced that to the order of (ka)? or (kb)*, ΠΠ = 0 and 





ο 


aa) Je = -£n 213 (ka)] {1-2ika ) "e ΠῚ 


n=] 


2ka 
= el = C + Lge oy 


MS etka 


x expl2 $ KV eu) fancı +» — AL 54] + 0(ka)? (5.15) 
gs V k*- wf 


By examination of the integral in (5.15) at the extremes of very small 


and very large values of u , it is shown that 


oo 


ku) {2n(1 + tk); du 


E Vk `= u 


is convergent for all values of k . Additionally, the integrable 
singularity at u = k presents no difficulty to evaluation. After 


substitution and a change of variable the integral becomes: 


fn mG x)9,(x) - AG (x)]+¥,(x)d (2 x) - (XY (2 de 
p [4,9 x) - 39.001 + 100-1007 
x {2n(1 + PENES dx (5.16) 
J1- (2) 


For small values of x integral (5.16) can be approximated by 
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ε 4.529 15.2 ER 2. 2 8 

In» Ee ut at, 
2 

0 (15 5) 


] b X b 
- dx £n x](2* 1) - 2n — 
z "ΕΕ --------------) dx 


b.e 
(4n 5) 
This last integral is finite, thus the integral (5.15) is convergent as 


u>0 . For very large values of x , (5.16) is proportional to 


oo A b ; 7b 3T T ] iti e 
| s sin(2 x - Fp) cos (x - 7) + 07 + [additional tems ]] 5 


Pie. x 
a SA 
k 4 
A L Fr) (2n(1 + a i 2 dx 
ο 4 ] 
| [1 αμ, (x) dx 


where we have used the expansion of ĉn(1+z) with z= A/x;z << l; 


z ž -l . Additionally the convergence of 


| sin mx ὦ dx and | cos mx (Š) dx 
M 


has been utilized. We are now able to conclude that the integral 


in (5.15) is convergent as u>œ, 
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Combining al] the foregoing, it is apparent that the integral 


appearing in (5.15) is convergent for all values of k and it is appar- 


ent from taking the limit as ka>0 of (5.16) that 


co 


lim | k Uu) {2n(1 ae), du = 0 


ka>0 k- u 


Additionally, the term involving the summation in (5.15) is very much 
smaller than one and therefore ignorable. The exponential factors are 
also negligible when large values of z are concerned. Therefore 


[x (k) T? in the low frequency limit becomes 
A = 2n u kae for kaykb << 1 (5.17) 
ao 


The reflected TEM mode field is now evaluated by substituting (5.17) 
into (5.14). The results are 


00 


He A = | Maly (ka) oe ik(zta)-iwt On (5.18) 
ποτ π΄, 


By the convolution theorem and wel] known integral (5.8.1) 


t: 
Ha (omz,t) = - | πι a) ale) dT zz (5.19) 
? ref lected TEM -οο 


where 





and 


roc lh 


If we again ignore for the moment the frequency restraints and allow a 


Er 
It} = 2 


Soe 


IA 

= 

IA 
oie 


Tr 


otherwise 


unit step function voltage excitation, then the reflected field is 


given by the integral of g(t). 


the following: 


me(p,zZ,t) = 
Po 
reflected TEM 
0 
€ 
- E == ) {5+ sin €t 
Ho Tp £n a e 2 
T 


and the electric field is 


= [5 Η t) 
E, pez st) EZ E p PZ? 


reflected TEM 


reflected TEM S 


In the case of a step input the reflected electric field is 


2a 
essentially the negative of the incident field for ut 1:1 J 


Carrying out the integration produces 


ol% 





Bone 
(see equation (2.6)). The current reflection coefficient is -J (ka) for 


ka << 1 , which is consistent with transmission line theory. For 


-t/T, 


allowable inputs of the form (1 - e with T large enough so that 


the low frequency restriction is not violated, the reflected 


ΠΤΙ 
fields are given by the above equations,multiplied by (1- e ra E 


ee a 
eet c Ξ 
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6. Radiation from the Pulse Excited Cylindrical Monopole Antenna 


The objective of this chapter is to examine the radiation from 
a pulse excited cylindrical monopole antenna in the light of results 
obtained from the preceding study of a coaxial waveguide antenna. We 
have determined that in the far zone the radiation appears to originate 
from the edges of the open ended structure located in the aperture 
plane. The edge starts radiating when the leading edge of the incident 
current pulse on the feed line impinges upon it. Radiation from the 
inner edge is opposite in sign to radiation from the outer edge. 

We begin by considering the approximate far zone fields from 
the annular slot antenna. This antenna has been the subject of several 
authors [28-32]. The structure of the annular slot antenna is that of 
a coaxial line terminated in an infinite nen plane as shown in 
Figure 6-1. Note that the center conductor is not hollow. 

We will follow the approach of Ref. [32] and replace the elec- 
tric field in the aperture by an equivalent magnetic current and 
develop the fields in the half space z > 0 from this equivalent 


source and image theory. The equivalent magnetic current 1s 


Since this source is just above a perfectly conducting ground plane, 


the electric vector potential is given by 








Figure 6-1. The annular slot antenna 
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l ik pe 
x € en b “Ε.Ε e | | 
F = = 0 p > i I i 
Jn | perc p dp dọ 
0 a Ir-r'| 
where 
nh > 
eE =-VxF 
(6 
and 
κ - τι - [o^* p'* - 200! cos($ - $') menr 


Now it is clear that the electric vector potential and the fields are 
independent of the coordinate 4 , from the cylindrical symmetry of 

the structure. Therefore we can set d& = 0 without loss of generality. 
We also know that the only nonzero component of F is the $ com- 
ponent due to the symmetry of the equivalent magnetic source. From 


geometry we find 


a, =a. cos ¢' - a. sin @' 
gy : x d 
since $70 ; ay = D for this formulation. Therefore 
27 b . 
e oL. p i a KR ! d "do! 
zm In ἃ 5 cos ϕ mp $ dp 
0 a 


with 


R= (p+ p'* - 2pp' cos ϕ' * (z-z')^) /* 


From Maxwell's equations and the Lorentz gauge it is an easy 


matter to show that 
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ο... E 32 
H = zu ΕΙ: s F] 


O 


Since F has only a $4 component which is independent of the angle 0$ 





> 
VeF = 0 
and 
iue. D en ΜΙ; 
ih - iu 2 = = >. | Es o'do' | Cos φ' en ἀφ' (6.1) 
a 0 


Equation (6.1) is exact--no approximations have been imposed. If En 
the electric field intensity in the aperture, were known exactly, then 
the integration could be carried out and the result would be exact. 

If we approximate the aperture field as being that of the inci- 


dent field in the aperture, i.e., 


Doou) ius 
o' gn = 
then 
: b 27 . 
uma ikR 
ge E "130 cos ϕ' "g— dp'do' (6.2) 
gn — 

a a 


In the far zone the distance R becomes 


R - (p^ * p'* - 2pp' cos à! * (z- z')^) /* 


- (r^s mes 2p'r sin6cos q' - 2rz'cos 0 * z'^) /* 


12 i 12 
r[1 + — - £D. sin 0 cos φ' - ££. cos 6 + ἔπ] 
r 


r 
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i 
wp] = on 0 cos ¢' | 
mance br o> o and r>> zZ, 


e KR B ur - p'sin 8cos φϕ') 


R B r 





We have retained more terms of the approximation to R in the phase 
factor than in the amplitude factor, since small variations in R 
cause large changes in the value of the exponent, The amplitude 
of the function is relatively insensitive to variations in R as 

r +o whereas the phase is not. 


Upon incorporation of the well known expansion [28, p.230] 


eikp'sindcos ¢ _ ) (-i] 2, (ko'sin SAULA 


= 00 


equation (6.2) becomes 





; : Ὁ ὅπ 
1we ikr oo NL. 
H = - τ Mole — | | cos 6' ) (-1)^9, (ko'sin ge do'do' 
f (Qn —)r nz-o 
a a 0 
: b 27 
ikr 
= V(o) e = | | 2 σος  ϕ' ἀφ' J, (kp'sin θ)άρ' 
u a 0 





£y 





we V(w) _ikr 2 d(J, (ko' sin 6)) 
| | k sin 6 

a 

r 








ik J (kb sin 8) - J_(ka sin 9) 
- fe EU ιν...” ] (6.3) 
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The approximate time behavior of the far zone field follows 


directly from the inverse Fourier transform of (6.3) 
m πι. 


e do . fe ESEL Üben Oi 
ϕ /21 Y Mo (gn 2)r sino 2 E 


- J {ka sin 0)]dw (6.4) 


When this integration is carried out by the convolution theorem, the 
impulse response of the annular slot antenna is found to be given by 
equation (5.10) and the step response is found to be given by equation 
6.11% 

From this result it is apparent that in the annular slot antenna 
tne radiation also appears to emanate from the edges in the aperture 
plane. The radiation starts when the incident current on the conduct- 
ing surface reaches the edge. One edge is located at p = b and the 
other at pop = a. For an incident impulse of voltage the outer edge 


contributes radiated magnetic intensity in the amount 


| NS 
[ 0 ΤΙ. ΠΠ | 
Z 
mes... Ee o -Psinec« «P2 sin 6 
V "o qr sin 62n > T er "a ac 
EC 
1 
; t> sine / 


(6.4.1) 
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and the inner edge contributes radiated magnetic intensity in the 


amount 
a 
0 ; T£- ¿SINO 
€ 
lb s ic - -$sinoxc «2 sine 
ο mr sin 6 £n = 2 sin ο. 
a 
0 ; T27 sin 0 
: r 
with t-t- c (6.4.2) 


We make the assertion at this point that the tip and the base of 
a cylindrical monopole antenna (Figure 6-2) give off radiation in a 
manner identical to that of the annular slot antenna. In this 
case, the radiation from the tip must be delayed in time by an amount 
which equals the time it takes light to travel the length of the 
antenna, h/c . We further assume that the currents on the conductors 
are reflected in the same manner as the currents on the conductors of 
the open ended coaxial waveguide. That is, the total current must go 
to zero at the tip and the base. 

In Appendix E we develop the radiated fields of a very long 

hollow cylindrical antenna. Transcribing those results, we have that 


the magnetic intensity of such an antenna is 


Eo Vlo) elkr J (ka sin 8) 
Ha tn 9,0) Em r sin 8 
o sin 8 £n = 





This result is identical to the frequency domain radiation from the 
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Figure 6-2. The cylindrical monopole antenna 
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tip of the center conductor of the annular slot antenna determined in 
(6.3) above. Since ka << 1, the particular form of the end does not 
have a significant effect on the fieldintensities. We therefore take 
these results to be valid for all thin cylindrical antennas regardless of 
the end configuration (with the obvious exception of those with end loadirg.) 
From the above discussion, it is clear that the contribution to 


the radiation from the tip of the center conductor is identical for 

h = Q and for extremely large values of h. We assume this result to be 
valid for all cyIndrical monopoles with h >0 . Since the base of the 
antenna is excited at t= 0 , when the incident excitation voltage 
reaches the ground plane, appropriate mathematical factors must be in- 
corporated into the field intensity equations to account for the 

finite time h/c required for the excitation current to reach the 

tip of the antenna. (c is the speed of light). 

Since the electric field at the junction of the feedline and 
the antenna must be continuous, when the flowing charges reflected 
from the tip reach the junction, charges of opposite sign must travel 
back along the antenna [8]. If we consider that the annular thickness 
(b-a) << A , then image theory is applicable and it is readily seen 
that the current relaunched onto the antenna from the base is equal 
in magnitude and opposite in sign to the current incident from the 
antenna. 

No radiation besides the initial radiation due to the current 
incident from the feedline originates from the base of the antenna. 
The current reflected from the antenna tip is cancelled by the image 


current at the z= 0 plane. 
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To determine the total contributions to the far zone radiation, 
use is made of Figure 6-3. Proper accounting must be made of the 
radiation from the tip of the antenna and that reflected from the ground 
plane. It is readily seen from this figure and the above discussion 


that the radiated magnetic intensity in the frequency domain is 


ο. ikr 
Ha (750,0) - [ο Viw) e [J (kb sin 0) - Jo (ka sin 0) 
Ho (2n Ir sin 0 


RIS - cos 0) , Q,IKh CI + cos 0) 
1 4 ezikh 


kh # (2n+1) 5 ; — n-0,1,2,3,.-- (6.6) 


The case for kh = (2n+1) 5 will be treated later. 
The time behavior follows from the inverse Fourier transform of 


(6.6) 


] -iwt 
(r,0,t) ες H (r,0,w) e dw (6.7) 
ó mm ) ὁ 


=O 


H 


If the excitation voltage is harmonic time dependent, then 
γ(ω) - βε[νΖπ A δίω- ως) 1 


From (6.6) and (6.7), we have 
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Figure 6-3. Effect of image plane 
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H¿(r,0,t) - [2 ——— Re IE sin 0) 
Ho (Xn it sin 0 


ia. cos 0) i oss 6) 
^s n8 ug og 


To t- L 
o iw ( e 


Upon simplification this becomes 


€- A cos w (t-f) w b 
y [So ONE ο 


o (£n zn sin 6 





ο C 
^os SOR 
C 


h 
[J (ο΄ sin 0) JLcos(— cos 0)] | 


woh = 
υπ ας ——4 (2n*l) » i ΠΞ0. 1.2.4... (6.8) 
ς 
This expression shows that the radiation pattern is 


wP woh Wa won 
A sin 8)cos Ser DE sin e) JLcos( — cos 6) | 
og Oe nn (6.9) 
sin 0 
The radiation pattern given by (6.9) can be compared with the pattem 
derived when the antenna is assumed to be very thin with a filamentary 


sinusoidal current distribution [33]. For that case 
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w h 


cos -τ- - cos (u ^ cos 0) 
F(8) * —————————————— (6.10) 
Sin 8 


In the limit when the feedline diameters are "thin" (i.e., A >> a,b)» 
the : patterns given by (6.9) and (6.10) are identical. Equation (6.8) 
shows that the radiation along 6=0 is always zero, as it must be. 
To get expressions for the cases where kh= (2n*1) m 
n=0,1,2,3,+..., we must consider the fact that ohmic losses are present 
in any real system allowing us to use a complex value of k or we can 


take k to be complex for the sake of mathematical expediency and let 


its imaginary part vanish at a later stage in the analysis. Thus, 


k = k, + ik k 


| “ο 
] 2 ? C 


ZEE 


With complex k the factor 


e ikh(1 - cos 6) " 2 kh (1 + cos 6) _ cos[kh cos 6 
Zıkh cos kh 


l+e 


cos [ (k,* ik,)h cos 6] 
cos(k,+ ik, 


cos(k.h cos 9)cosh(k,h cos 0) - i sin(k,h cos 6)sinh(k,h cos6) 


COS Kıh cos ADS sin 1 sinhiK, 


In the limit as kh > (2n+1) = n=-0,1,2,3,°**, this expression is 
seen to be finite for Κρ #0 and the magnitude is only limited by 


the fact that the conductivity of the material is not infinite 


= ED" d | 
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TT 
cos(kh cos 8) cos [(2n*l)s cos 0] 
cos k 


-iLsin(2n*1)5]s inh (k h) 


lim E 
kih > (2n+1)5 2 


n=0,1,2,3,--- and Κρ -ε»0 


The magnetic intensity becomes 


B Α οἰπίω 1 - Ὁ 
H¿(r,0,t) ~ [ο ------------ 
V Ho (£n μπα 


ao. T 
ο ος sin 9)[cos((2n+1)7 cos 8) ] 


sin[(2n*1)5] sinh(k,h) 


The input current to the antenna peaks at the values of kh 
given by kh = (2n+1)> , n=0,1,2,3,°°° . The explanation for this is 
that the monopole can be considered as a simple extension of the feed- 
line. The tip of the antenna would correspond to the open end of a 
transmission line. Any odd number of quarter wavelengthsfrom the 
tip along the line toward the source would correspond to a point where 
the impedance is a minimum (a point of series resonance) and the cur- 
rent is a maximum. When kh = (2n*1)5 , the feed point is an odd 
number of quarter wavelengths from the open end of the extended trans- 
mission line. 

Taking the inverse Fourier transform of (6.6), we find the 


radiated magnetic intensity for a voltage impulse of excitation 





-110- 
] 


. [5 |(ς sin 0)^- 1 
Ha (7,0, t) - i 


: 0«c«2 sin 6 
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impulse man T)r sin 6 0 ; elsewhere 
exci tation 
Tora Ξ 
C 
co n 
( ) νει -- >, O0<2 <2 sine 
120 J (X. sin e)? e ane: 
P C N 
i 
| 0 ; elsewhere 
E rs hy, | . enh 
πη Cu ο cos 6) ee 
Bi n 
EM . ues 
no uno 22 
f ' (= sin 6) T 
0 s elsewhere 


ce Eh nh 
DU ELE lit COS 6) - Ξ (6.12) 


For a unit step of voltage excitation the radiated magnetic 


intensity is the integral of the preceding expression. 





ie 


H¿(r,0,t) 


unit step 
excitation 
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C 
| 
€ 
Ae — H+ sin I); - Bsin O<T<T SINO 
Fo lan =)r sin 6 \ 
| 
| T T > D sin 8 
C σος 
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" r h 2nh 
f ee ul cos 9) Ξ 
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= TIC 
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(6.13) 
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Excitation Voltage (c) 


(volts ) 
Loc 


(a) | — (b) 
ya 









b sin 6 





b sin 8 
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\ / b sin : 


Figure 6-6. Radiated magnetic intensity of annular slot antenna 


i b. 
Ordinate: Ha = ΝΡ ΓΗ, cat 9)(volts) 


Antenna Dimensions: h= 0 and a=. &b 
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The magnetic intensity for arbitrary time dependence can be ob- 
tained by convolution of the input excitation voltage with the impulse 
response (6.12) or by convolution of the time derivative of the input 
voltage with the step response (6.13). As expected, the radiated field 
for the infinite monopole is the same as that for the finite length 
monopole for the time interval before the current pulse reaches the tip 
of the antenna. The results, equations (6.12) and (6.13), are easily 
extended to the infinite cylindrical antenna by letting h >œ . The 
radiated fields for the annular slot are found by letting h » 0 in 
expression (6.6). It is also worth observing that the results (6.12) 
and (6.13) are based on the assumption of no ohmic losses as the cur- 
rent pulse travels back and forth between base and tip of the antenna. 
The attenuation can be found using the conventional transmission line 
approach for small attenuations on good conductors. The determination 
of the time domain fields in the case of a lossy conductor is com- 
plicated since the resistance per unit length is proportional to the 
Square root of the frequency. The incident pulse supplies the total 
of the radiated energy, ohmic losses, reflected energy and energy 
stored in the jnduction fields. The stored energy of the induction 
fields is the source of energy for the second and subsequent pulses 
launched down the feedline toward the source, the energy radiated after 
time 2h/c, and the ohmic losses which occur after time 2h/c . 

The theoretical results obtained in this chapter are in good 
agreement with the experimental and numerical work of Schmitt, et al 
[6], Palciauskas and Beam [10], King and Harrison [8], Burrell [34], 


and Lamensdorf [14]. 





-116- 


Figure 6-7 is a comparison of the theoretical results obtained in 
this work with the experimental results of Schmitt, et al. [6]. Reference 
[6] states that the second pulse measured in the experiment should be 
moved ./4 ns toward the origin and its amplitude should be about 3/2 as 
much as observed. This is necessary to account for experimental devia- 
tions associated with locating the receiving probe a finite distance from 
the antenna. In Figure 6-7 these corrections are applied to the second 
and subsequent pulses. 

It is also constructive to include the following observation. 
When a cylindrical dipole antenna is infinitely thin and center fed, it 
is common to consider it as an open ended transmission line that has 
been spread or opened out. A sinusoidal antenna current distribution 
with current nodes at the ends is assumed and the radiated fields are 
obtained on that basis. If it is further assumed that the current at 
the antenna feed point is equal to the current existing at this point 
on an equivalent length of open ended transmission line, then the field 
expressions developed from the assumed sinusoidal current distribution 
are identical to the field expressions obtained in this work in the 
limit of very thin antennas. 

To demonstrate this point, we first observe that the radiated 
magnetic intensity created by a filamentary sinusoidal current distribu- 
tion of I sin k(h - |z|) is [33] 

RT 


[De [eos kh - cos(kh cos 9)] 
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The driving point current is I sin kh . The other symbols have the 
uSual meanings consistent with their previous use within this chapter. 
It is also common knowledge that the input impedance of an 


open ended lossless transmission line of length h is 


Z = -1Z cot kh 


where Lo is the characteristic impedance of the line. If we drive 
such a line by a perfect voltage source V(w), the source current is 
readily obtained by Ohm's law. Since the antenna driving point cur- 
rent is assumed to be equal to the current which would exist at the 
input of an equivalent length of open ended transmission line, we 
equate the driving point current to the source current. Therefore, 
I sin kh = πας 


and 
V(w) 


lo Γ -12 cos kh 
By substituting this value of - into the field intensity equation, 


we derive the general result for very thin antennas, 


_ 2-V(u) elkr ΠΕ cos(kh cos 9) > 


Hr, 9,6) - 2π7 Y sin 9 cos kh (6.14) 


Outside of a discrepancy in sign, it is an easy matter to show that 
this general result reduces to equation (6.6) when the limiting form 
for small arguments replaces the Bessel functions of equation (6.6) 


and the characteristic impedance of a lossless coaxial line is 
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substituted, 1.e., 


om 
iI 
INS] — 
3| 

01 E 

o [o 
2 
= 

£& |c 


The discrepancy in sign is the result of not specifying a 
voltage reference point and it is not significant relative to the 
results. 

The approximate transient behavior of a very thin dipole 
antenna fed by a transmission line with characteristic impedance 
Z is found by taking the inverse Fourier transform of equation 


0 
(6.14). 
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7. Conclusions 


The results presented in this report provide a comprehensive 
analysis of the open ended coaxial transmission line or coaxial wave- 
guide antenna. The fields radiated from such structures can be com- 
puted without a priori assumption of the aperture field. The current 
reflection coefficient for the fundamental mode is found to be 
-J (ka) for ka << 1 

The conventional approximation for the fields radiated from 
the coaxial aperture antenna gives good agreement with the results ob- 
tained through the Wiener-Hopf method for thin radii in the feedline 
(kb,ka << 1). 

When the excitation is such that A << b-a, the radiation 
pattern obtained by the conventional analysis method is shown to be 
valid for the predominantly forward direction only. We determine that 
the classical radiation pattern obtained by conventional methods should 
be multiplied by the factor [cos | 1/2 for 1>0. The resulting 
pattern is one in which the radiation at high frequencies is confined 
to the forward direction as expected from geometrical optics predic- 
tions. 

A model for the cylindrical monopole or dipole antenna is 
developed based on the observation that radiation appears to emanate 
only from the base and tip of the structure in pulse excitation. This 
model is used to determine the radiated electromagnetic field of 
cylindrical monopole antennas for both harmonic time dependent and 


arbitrary time dependent excitation voltages. 
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The mathematical techniques required to obtain results using 
this model are EN simple and in most cases hand calculations 
are sufficient. The results obtained show good agreement with experi- 
ment and are useful for both analysis of the radiation field when the 
driving voltage is given and for synthesis of a driving voltage to 


produce a specific radiation field. 


Future research can be aimed at removing the restriction that 
the width of the annular region of the coaxial waveguide antenna is 
very much smaller than the wavelength of the exciting source. The 
results could also be extended to the thick antenna case (i.e., 
kb,ka » 1). Additional research effort on the transient behavior of 
the cylindrical monopole could be directed toward characterizing the 
effect produced by different end cap confiqurations, determination of 
frequency limitations, accounting for dispersion and describing the 
effect of hiaher order mode components of the aperture field on 


the transient radiation. 


ο... 
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Appendix A - Analytical Determination that H Modes are not Present 


on Structure 


In this appendix we will show that H modes are not present on 
the coaxial waveguide antenna excited by a TEM mode. 
It is obvious that symmetry precludes excitation of all H modes except 
ο (n=1,2,3,°°:) modes. 

By the edge conditions [19] (see Figure 2-1) 
we know that the longitudinal current must approach zero. Thus the 
sum of the incident current and the reflected (longitudinal) current 
must be zero. Therefore the reflected current on each conductor has a 
longitudinal component equal and opposite to the incident current. All 
of the reflected current is longitudinal as it must be to completely can- 
ce] the incident current. To have jos modes, a circumferential com- 
ponent of current should be present on the structure. Since no current 
source is available to launch a circumferential current, the Um modes 
cannot be present. 

The above paragraph notwithstanding, let us assume that some 
current is launched in the circumferential direction upon reflection 
of the incident current. It is shown in Chapter 2 that the longitudi- 
nal current in each conductor is equal and opposite. In keeping with 
this, we assume that the reflected currents on inner and outer conduc- 
tors must also be equal in magnitude and directed in opposite circum- 
ferential directions. It will now be shown that me 0 (n:1,2,3,:--). 


The applicable scalar Helmholtz equation for this problem is 
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yt tye lee = Eg (052,0) = 0 | (A.1) 
Es (02,0) = 0 at p --^asb sz «9 
Eq (092,0) = 0 at p=0 3; <z<o 
with edge condition: 
+ 


SAID +0 at p.=ab; z>0 


Fourier transforming (A.1) gives: 


2 


y t -γ΄- 7 E (Pasu) - 0 (A.2) 


D |= 
vja 
èl 


Using the regions of Figure 2-2, with H, replaced by E, and 


$ $ 
the asymptotic behavior of E „(p»2,W) as z+ +o , we find that 


Ε, (ρνανω), E, (p,a,w) and E, (p,a,w) are analytic in the strip 
$1 $2 $3 
shown in Figure 2-4. Solving equation (A.2) with Es = Es , Es , Or 
1 2 
E, respectively, yields 


$3 
ΓΦ) Ξ AK4 (vo) (A.3) 
Eo = BI} (yọ) + CK- (Ypo) (A.4) 
7 = DI, (yo) (A.5) 


(A.3) satisfies the radiation condition and (A.5) goes to zero at 


p=0. At p=a, 
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Εφ (άνανω) Ξ E, (sou) (A.6) 
nen) = Eg, aseo) Ξ0 (A.7) 
d T d + 
EEE “Ν σος. (o,0,u) ] B (A.8) 


Eg, (rca) = BI (ya) + CK, (ya) 


or 


Ερ(άνα,ω) - Β{1η(γ8)) 








ς = ΕΠΟΝ (Α.9) 
d Toe 1 B Et (oa) A. 10 
dot? g, (Pau - = Κίνα] ` ya πιο (Α. 10) 
From (A.5) we have 
EY (a,a,w) I (ya) 
-[οΕ (ρ»α»ω) ] Ξ γᾶ "3 i (A.11) 
p 9 Oy = ae . 
p $4 o=a I; ya 
where use was made of 
E (a,a,w) 
_ 93 
D Πτα (A.11.1) 


Subtracting (A.11) from (A.10) and using (A.8), it follows that 
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d - d _ 
ae, (ρνανω)] nd a (Pu) 


p-a 


Kg(va)  Ij(va) 
a rt rn 


Define the left hand side of the above equation.as G (a) 


E? (a,0,0) 
- $ 
G (a) = B + 3 
K, (ya) ΚΙ γᾶ 1 ya 


At p=b: 
+ A 
Ey, (bau) Eg, (bsc) 
u - Egy Dect zo 


d + d + 
pE (p,a,w)] = ~ pE (p,a,w)] 
d$ pb ορ 42 p-b 


Eg, (Peau) = BI (yb) + CK- (yb) 


2 Bella) = BI, (vb) 
9 
C z 
K, (yb) 


Str, (ao) x = yb[BI (yb) - CK (yb)] 


B K (yb) 
"KG5)- - yDE, ἕνα; ω) K 5 


(A.12) 


(A.13) 


(A.14) 


(A.15) 


(A.16) 


li) 





s ] το πο ο A.18 
dpt? e, Pi) m m 9, 9%, W Ky (yb) ( . ) 
where we have used 
(b,a,w) 
E ] 
A = πα (Α. 19) 


Subtracting (A.17) from (A.18) and using (A.15) it follows that 


d = d > B 
pE (p,a,w)] x [oE (0,0,0) ] um 
d^. ob «Ὁ $ b Kb 
Define the left hand side of the above as E (a) . Then 


EHE KUBY (A.20) 


If we substitute (A.20) into (A.12) we get 


+ 
E Κι (yb) Eg aan) 
(9) - Ce) ray * Ino — 


By equating the values of C given in (A.16) and (A.9) , it follows that 


Ar K, (ya) I (yb)  1,(ya) 
iau) 7 Kay Εφ (αμ ανω) - BL TET Πε] κ K, (va) 


Upon use of equation (A.20) the above relation reduces to 


+ K, (ya) & _ 
Eg s) ΙΤ Eg, sao) = E (a)[K, (ya)1, (yb) - 1, (ya)K, (yb)]. 


Substitution of this equation into (A.21) gives the companion equation: 
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E* (b,o,u) 
G”(a) - E (a) m" es (A.22) 


From the definition of G (a) we note that it is proportional to 
H_ (a,a,w) - H_ (a,a,w), and E (a) is proportional to H_ (b,a,w) 
23 22 2] 
EM (b,a,w) with the same proportionality factors. By well known 
2 
boundary conditions we know the surface currents on the outer conduc- 


tors are, respectively: 


> + E 
~ - = + 
a XL, (20,0) N, (aw) Ja, Ly a, 


i 
J 
w+ 


a AER, (bos - κ... Y 


In this problem we have assumed that b-a << ae and therefore the 


parallel plate case 1S approached. Under these conditions with TEM 
excitation, the currents on the plates can be taken as equal in mag- 
nitude and opposite in direction, with a great degree of reliability. 
Therefore G (a) = -E (a) and equation (A.21) becomes 


Es (a,a,w) 
: 3 
G (a) * —————————————— (A.23) 


I (va) [K, (yb) + K, (ya) ] 


Simultaneously (A.22) becomes 


E (b,a,w) 
A AS (A.24) 


K (yb )LI, (yb) + I (va)] 


We may factor the denominator of the right hand member in (A.23) and 
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(A.24), using the techniques of Appendix D. 


K (yb) (1, (yb) + J,(ya)] = N (o)N (o) 


with Na) = N"(-a) « (lal V^ as aro 
with |Im a] » “ko 
Simi larly 


1, (ya)[K, (yb) + K,(ya)] = M(a) M (a) 
wi th M (a) = M (-a) a ΠΠ. as ato 


with [Im a| > -k, 


It is also necessary to point out that G (a) % a. as a9, 


This is determinable from the edge condition that H,(0,Z,w) v iz 
dez > 0 at p=ab. 
Substitution of these factored forms into (A.24) and (A.23) 


and algebraic operation yield the entire or integral functions (A.25) 


and (A.26) 
E, (b,a,u) 
G (a) M (a) - 7 (A.25) 
α 
E (a,a,w) 
G (a) N (a) = — (A.26) 
α 


Each side of (A.25) or (A.26) is the analytic continuation of the 
other side. Both sides of (A.25) and (A.26) are analytic in the strip 
-k5 ET « k, . Since each side of (A.25) and (A.26) is bounded for 
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all a, they must also be constant (Liouville theorem). Examination 
of the asymptotic behavior of the above expressions shows that each 
side is identically zero. ~ (a,a,w) +0 as lal »» , since 
3 

E.(a,z,w) >0 as z> 0°. Therefore El (b,a,w) = X (a,a,w) = 0. 
$ >) $3 
Each of the coefficients A,B,C,D can be expressed by a linear combina- 
tion of E, (b,a,w) and S (a,a,w). This conclusion is evident by 

1 3 


observation of equations (A.11.1), (A.19), and the simultaneous solu- 


tion of 
EN (a,a,w) = BI,(ya) + CK,(ya) 
$3 1 1 
E, (b,a,w) = B1,(yb) + CK, (yb) 
by 1 1 
Hence E, ,E and E, are identically zero and E,(0,2,0) does not 
9° ὁ) Φα $ 


appear in any region of the coaxial antenna structure. 
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Appendix B. Determination of Branch Cuts for y 


In this appendix we will develop the branch cuts in the com- 
plex o plane that define the region of analyticity of Y , a multi- 


valued function of a 
y =/a-k 


Our starting point for the work in this appendix will be equations 


(2.9) and (2.9.1): 


2 


d 1d ] 2 
A - - vy ) H, (ρνανω) - 0 (2.9) 
do p do o? $1 
H, (o,0,0) » 0 as pre (2.9.1) 
97 


The solution to equation (2.9) satisfying (2.9.1) was found in Chapter 
2, 


He (p,a,w) = Ala) K,(ro) (2.15) 
1 


The asymptotic value of the modified Bessel function [20] is known to 


be 


E CD 00 
K- (Yp) l e as p> 


Hence 





He (ρ»Ζνω) œ l | Ala) e YP e719% da 
Mel 


as opo (B.1) 
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From the Sommerfeld radiation condition we know that the field inten- 
sity must be an outward traveling or evanescent wave for large 
r= jot z 
Following the manner of Reference [21, p.20], we note that if 


we require that 


Re y > 0 (B.2) 
Imy <0 (B.3) 


then the radiation condition will be satisfied by Hy (032,0) . This 
may be verified by examination of (B.1). The branch cuts will be 
chosen in a manner which ensures that the conditions (B.2) and (B.3) 
are satisfied. 

The a plane is viewed as a two-sheeted Riemann surface with 
the sheets connected along the branch cut. In each of the sheets the 
function y is single valued. The sign of the function on one sheet 
is the negative of the sign on the other sheet. The value of y 
becomes discontinuous only if a branch cut in the a plane is crossed. 

Let us define the branch cut of y such that Rey >0O in 
the top sheet and Re y <0 in the bottom sheet. Thus the two sheets 
are joined together by the curve given by Rey=0. 

The branch point is obviously k and we will, for computa- 


tional purposes, take k complex 


k= κι + ik, 3 S >> Κο : Kısk, > 0 


We may therefore write 
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2 pu x 


EL = (c + it)* - (k, + ik) 


(0° - 3^) - (t - KÀ) 4 21(στ - Κηκρ) (B.4) 


Note that Re y > 0 in the entire top sheet, therefore if we 
specify E - re! and y = p 1/2 10/2 with 0 <8< 2r , then 
Re y > 0 only if lo] «vr . Thus the branch cut must be given by 
|o| = T or Re 2:0 and Im r = 0 . The branch cuts which divide 
the a plane can be easily identified by examination of Figure B-1 
and are seen to be a portion of the curve oT= kk, : 


Since the branch cut is given by d = re!" it follows that 


1/2 ,im/2 or Rey 2-0. The graph (B-2) depicts 


on the branch cut y=r 
the branch cuts for the top sheet. 

For the branch cuts in the bottom sheet, it is only necessary 
to use the opposite sign of y everywhere in Figure B-2. 


It is also worthwhile to state that for this branch definition 


$17 6 
y = D Εν where D —— 


A p ad 2947 
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zl E 
2, Im y` = CT i. 
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Figure B-1. Regions of complex o plane 
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Figure B-2. Branch cut for multivalued function y 
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Appendix C. Useful Expressions and Expansions 
In the main body of this work,extensive use is made of ordinary 
and modified Bessel functions and Hankel functions. Listed below are 
explicit expressions, asymptotic approximations and limiting forms for 


some of these functions [20,24] of integer order. 


(i) Relations between modified and ordinary Bessel functions: 


1 : pe 
- πνπὶ m 
e 2 J e? Ίο (-π < arg z en 


I (z) 


Mj — 


1 à ]. 
vri mi 
ni e? nC ae? ) 


K (2) y (=m<argz<>) 


ES = J (z) + j Y, (2) 


(2 E : 
H es = y, (z)-1Y,(z) 


(11) Wronskians: 


NÍ— 


το τε ο ES 


Li l2) Y (2) - δν (2) Y, 7 E 


(111) Series: 
rz (4 


J A 
of) (11) (2!) 


ye (4 3 


τ ... 
(31) 





Nc Ic? 
ο παπα. 


(iv) 


(v) 


Limiting forms for small arguments (when v is fixed and z +0): 
J (z) ~ ίνα) ον. νι. 

¥ (2) N -Hf 1) (2) N in’) (2) N £ {η 7 

ya iD a infz) a - MY) (FY; πον) » ο 


M (> 2)"/T(v+1) ;  v # -1,-2,-3,--- 


K(z) ~ -2n z 

] ζΖι-ν 
K (z) 2 τν) (2) 
Asymptotic Expansions (when v is fixed and |z| >œ) 


J (z) a= [cos(z - 7 VI - 7)) € Ap αρα εἰ « T) 


Y (z) = JÆ Isintz - 1s - 5) « oip la ee 


Z 
NE (|arg z| « 2) 
ν JITZ 2 


TZ 





ar fetos (arg al <p) 
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(vi) Derivatives: 
d r_v ZN 
jn ολ... Luz) 


Sz KG» -2"K, i) 


(vii) Recurrence Relations: 
Ia ιν. (z) = 21 (2) 
ν-] v+] ZEN 


2V 


Kyo 1 (2) - K 41 6) Ν zuge 


(viii) Other Relations: 


I (z) 


=v 


D (z) 


K 2) - K (z) 
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Appendix D. Factorization Procedure 


The purpose of this appendix is to furnish a factorization pro- 
cedure for the Wiener-Hopf kernels which arise throughout this work. 
The factorization procedure used extensively in this work and developed 
in this appendix was originated by Bates and Mittra [23]. The mathe- 
matical development of this procedure is first presented and then it is 
applied to the specific expressions to be factored. 

We begin this appendix with a factorization of me wi th 


y = "a 2 using the well known theorem of [4]. 


Theorem A: Let f(a) be an analytic function of «a = gtit regular 
in the strip 1_<T<T, , such that fla) < clol” ,P>0 for 
|o| + , the inequality holding uniformly for all qt in the strip. 


Then f(a) can be decomposed such that 


f(a) = f (a) + f (a) 


with 
œ+] C 
f la) = X f(8) dg : NES ο κ τς τ 
rc pa 
oct i d | 
= E ] {{6 : 
f (a) = - 5 | FB) ag ; ET «ed s. 
-œ+ } d 


where f'(a) is analytic in the portion of the a plane defined by 
τὸτ and f (a) is analytic in that portion of the a plane de- 


fined by T< T}. 
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Theorem B: If £n g(a) satisfies the conditions of Theorem A, im- 
plying that g(a) is an analytic function of a , which is regular 


and nonzero in the strip 1 <zT<T, and gla) > 1 uniformly as 


+ 


o + +o in the strip, then g(a) can be factored such that 


g(a) = g (a) g™(a) 


with 
octic 
g (a) = exp | 218) qs o CK T. 
-+| C 
ovt id 
g (a) =exp - an | in ale) gg A SA 
-tid 


The function g (a) is analytic in that portion of the o plane 
defined by t >t and g(a) is analytic in that portion of the 


plane defined by t< tT The principal branch of the logarithmic 


D 
function is used throughout. It is also worthy to note that if 


ala) a ar and £n g(a) v» -Anlal as Ja] »» , Theorem B can still 


be applied. The integrals are convergent in the sense 


T+it 
lim | { } dB 
> -T+iT 


If f(a) and g(a) are even functions of a then 


fl-a) = fla) or fla) = fl-a) 


and 





E 
g’(-a) = g(a) or  g(a)= g(-a) 


The function e Yh , which we wish to factor is interpreted as 


e^ Ja - k co at » k* au JKT- a 2 2 | 


or for k >a Examine 


an g(a) = -yh = -h(o*- k*)(a*- k^): /* 


The first factor is an entire or integral function of a . Therefore 
we only need to decompose the second factor 1/y . Since y is mul- 
tivalued, we will use the branch cuts for the top sheet of the two 
sheeted Riemann surface of the a plane (see Appendix B). 


From Theorem A, we have 


oic 


f (a) = gr A 
3 P ht- k` (B-a) 


The integration path for equation (D.1) is diagrammed in Figure D-1. 


: -K, ee (D.1) 


The integral (D.1) will be evaluated by integration along a contour in 
the upper half plane, which does not cross the branch cut located in 
that half plane (see Figure D-2). By Cauchy's theorem the value of 
the integral (D.1) is 


fla) lo wee | — (D.2) 


dg 
Q = k Pai P, vB - k (B-a) 


The contributions to the integral from the semicircular portions of 


the path at .|8|» » are zero by Jordan's lema. 
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ı Imß 
B plane 
iko k 
N 
pe 
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N iR 
-k 
k = k, + ik, 


Figure D-1. 
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Figure D-2. 





Figure D-3. 
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For the limiting case of no losses, κ, 2 ο”, the path nen P, 
is adjacent to the coordinates as shown in Figure D-3. Under these 
circumstances the remaining integral of (D.2) is expressed as the 


five separate integrals which follow: 


0 k 
i — 


2 dx 
i Vya k* (iy-a) | i Jê- x (x-a) 


=T ir a! 0 do 0 dx 
+ lim go NF N à | 
ΤΉ ΚΓ ΠΕΙ ΚΠ ἘΝ 


a i dy 


-i Vy + k^ (iy-a) 


By reversing the limits on the last two integrals and evaluation of 


the third, we can reduce the above set of integrals to: 


k œ 
gali a , —— 
T" T 
0 NL. x 0 dy? Ke (y+ia) 


A further reduction is achieved by a change of variable. Let 


=ksins and y= ik sin 6 . Incorporation of this change yields: 


7/2 -jœ 
T l1 | —— 
7 0 k sin 6 -a y 0 k sin 6 -a 


Note that use has been made of the identity, sin(-iy) = -i sinh y . 
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From Reference [20, p.78] we have the following integral: 





D 
a tan 5 + b 2 2 | 
SATZ Y παν 77 ] NM 3 - 
à 
| dz _ 
atb sin z — Z +b- (bé. ac) V? 


1 2 7 
37935 10 > a”) 
(b A a tan 5 *b *(b e ac) i 


Through this integral we are able to complete the evaluation as fol- 


lows: 


- Q&Q tan 8 4.67 -atan 7+k 7 
arctan αι απ ae : 
Fm ee E 37 | (al. τα] i 


-Q tan 5 ας 2% 2 = 5 


α tan z*k-( B 
BECA i - &n 
τσ μα a aF u. — 5+ kt(k d yt 


0 


(ké > al) 


Since tan(-iy) = -i tanhy and lim tanhy = ] , the above relation 


y >0 
simplifies to: 
΄ 
2 en] 2 2 
arctan ( ) ao) 
NE | 3 
} 
| 1 *ia*k -(k dn -at k (κέ- 2 
aa ee 
5 
n (k^ -o^) / E *t (k^-a y? EET yu 


1% > 
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The trigonometric identity 


+ 
arctan z, + arctan z, = arcan 2] 
1 + 2,2) 


was applied to the top expression. 


Recognizing the relation [25, p.62] 


2d. l - iz 
> oe 


and applying it, we are able to achieve great simplification 


1 - (atk 2 2 
XY VEER s Gk 


Y = (at k)> o*) 


Em 
ty "hey > k> 


2 2 


Note that the same result is found for αἵ >k and k” > - . When 
the above is reinserted into (D.2), the decomposition is completed 


and 


K 


Pr + Sn 

2 emi > REE 

TY Qo toO 

tU PUE © ay UOTE 


The factorization of g(a) = ¿YN can now be easily completed. From 


a) 


the above we have 





2149 - 


ihy 


en g (a) um ee nfl + (a-k ] 


= QU -~ 


Because g(a) is even in a , it follows that 


+ 2 
n g (-a) = £n g (a) 
and 
^ ., ih - (a*k 
an do) DR n [ΠΗ 
Finally 
+ iby 


g(a) = g(a) =e "o mt (0.3) 


This result is consistent with that of Noble [4, pp.20-21], and is 


shown to be correct by the following: 


ΥΣ - (atk); 
g(a) = g (a) g (a) =e a E: nly t (am 


= € = e^ hY (D.3.1) 


where in z= Anlz| + ig 


-M< QST 


The only possible singularity of g (a) in the upper half plane is a 
branch point singularity. It is obvious that this singularity is not 
present and that g (a) is regular in the upper half plane, since 


substitution of si|y| for y has the same result. Similarly, for 


g (a) , 
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Y ll 1 
ein age. +ih fen e 197 { -h$ ζ 
e An[- ]= e -e 
| | Y ON 


- τα, 
9 = 2 tan oy 


Theorem C: [23] Let G(a) be an analytic function of a in the 


region tT <a<t, with the following conditions: 
1. G(o) is regular in the region tT <a<t,. 
2. Gla) #03; Gla) = G(-a) (even). 


ν -Π[α] 


3. Gla) v ca” e as |Reao| » where v and h are 


real constants. 


Then G(a) can be represented within the strip by 


G(a) » G'(a) G (o) 


where G (a) is analytic and nonzero in the half plane t»: and 
G (a) is analytic and nonzero in the half plane t<t,. The 


expressions for Bo) and G (a) are given by 
G (o) = G (-a) = { /G(0) (1+ v) 


(8* (a) exprime 4η ren ; ikh] } (D.4) 


with 
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(a) οο- ὦ 

+ | F 

B (a) = oo| zir | s sE) s| >; [Im al < [dl]; |t_| > ld] 

S (D.4.1) 

(b) 

F(8) « h(g^- k^) V? - v sn(g^- k^) V^ - on c « en cg) (D.4.2) 
(c) 

y = co UL with branch of Appendix B. 
(d) 

k = κι’ ik, 3 Κ᾿ >> Κρ} Kook, uw |r| <k, I. < k, 


We will now show the validity of the above theorem. We begin by observ- 


ing that 


: y PRO -aF(B) 
B (-α) = exp Εἰ | Bata) dg 


~o-j d 


If 8 is replaced by -B and use is made of the fact that 
F(ß) = F(-8) then B (-a) becomes: 


-+j d 
] αξ(β) 4β 
- È | sas | 


eot i d 


Β΄ (-α) 


and 


a βίο 


We next find the product 
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B*(a) B(a) = exp p: | te) | 
ς 
where c is the contour shown in Figure D-4. Observe that the asymp- 
totic behavior of the integrand can be used to show that there is no 
contribution to the integral from the vertical segments of the contour 
at |Reg| >» . By Cauchy's theorem, the value of the integral 


exponent in the above expression is 


21i ) (residues of the integrand, from the poles at 
B=0 and B=a) 


-F(O) + F(a) 


hy + ihk * £n | 
6(0)(1 2) 


Therefore 
B'(a) B^ (a) » Ga) exp[hy = 
G(0)[1 - =] 
k 
From (D.4), we formulate 


ν/2 


+ - E t - 
G (a) G (a) = G(0)(1 - 2) B (a) B (a) 


exp | ant, i sa-k ] + BY np - η 


Substituting for EN B (a) and using equation (D.3.1), we find 


" 





-15]- 


i Img! 


Figure D-4. 





ων" 
| 
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G (a) G (a) » G(a) , which establishes Theorem C. 

Precise knowledge of the asymptotic behavior of G (o) = G (-a) 
as |o| » within c» 1 is necessary for successful use of the 
Wiener-Hopf method. As our next task, we will determine this asymptotic 
behavior. Notice that μμ. is bounded as lal +œ , therefore by 


letting [al +o in (D.4), we find that 





‚a + (a-k) 
1 Qn 
EUN SEC UE Clay) line MEC. 
Q | +00 
. ah δα 
τ. 
for |a| » within 1 > _. In the manner of [23] the factoriza- 


tion formula of Theorem C will now be converted to a form convenient 
for numerical work. In this development, we concentrate on changing 


+o-1d 


+ _ ] aF(8 
B (a) = exp E | B(B-o e 
-o-id 
to a more tractable expression. 


Taking the logarithm and integrating by parts yields 


οο-.1 


2n B (a) = Jy ja F(8) (ον) € 
de: 

= F(p) an(B3)-[ Α(Β) and) ag 
00 d 


where 
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A(g) = Bh ο βν 2 ] dG(8) 
Zm(ke-ge) 7 οπι(κ2-α2) mi ce) E 


If we let B> -B , then 
—o+id 
in Ba) = | AC) sn(72 ag 
cot] d 
This integral will be evaluated on a closed contour in the upper half 


plane by the residue theorem. Therefore 


en B (a) = - | A(-8) &n(1 + 5 d& - 2mi ) (residues at the poles 
of integrand 1nside 
contour ) 
The path p and the contour are given in Figure D-5. The contribu- 
tion from the semicircular path at |8| >œ is zero. 
We now assume that 
(1) G(a) has at most one branch singularity in the upper half 
plane located at a= +k in the fom y= lal- key Ve. 
(2) G(a) has only simple poles in the upper half of the a 
plane, located at a = Τρ with | Im Pr, | > T, 5 n=1,2,3,--- 
(3) G(a) has only simple zeros in the upper half of the a 
plane located at a = z, with | Im z IS, 
A TNI of A(ß) shows that its only possible poles within 
the contour are those of ata SS The poles of this function 


occur at the poles and the zeros of G(8) . The residue contribution 
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ea Reß 


Figure D-5. 
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from the poles of G(8) is GE) and from the zeros of G(8) it 
is (gr). Since G(8) is an even function, when 8 is replaced by 
(-8) the contributions from the poles and zeros become the negative 


of the above. We may now write 


gn B (a) = - | ASI A A ES 
n=] Pn πε] Zm 
p 
Set 
a ETE 


where the branch on the upper Riemann sheet is used (see Appendix B) 


β - i lu- k2) V2 _ (κέ- ut 1/2 
Using this substitution, we have 


4 
μυ --- | A(-8) an(1 + 2)de = KD anf ume 


p p κ. 
where 
Cu) = + gv - XT E [2n G((k?-u£) 23 (D.6) 


(see Figure D-6). 

From Appendix B we know that on the branch cut (g*- κε) ED. 
therefore (gé- k*) is pure imaginary and the variable u is pure 
real. So the introduction of the variable change maps the contour p 
into a contour along the real axis. The contour p' must be in- 
dented around any poles of the integrand which appear on the real 


axis. The indentations are consistent with pole movement from the 
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real axis if the medium is assumed slightly lossy. Note also that we 
have accounted for a negative sign by taking p' in the opposite 
direction to the path which p maps into under the variable change. 
Examining C(u), we note that the second term has a pole at 
the origin,and the third term 


-ilda ο. ο οι... 1 d 2 2,1/2 
Ani du [an G((k -u )] (ami) (κ -γ du G((k^-u ) 


could possibly have a pole at the origin, and at other locations along 
the real axis, which follows from the fact that u is real. The 
contribution to £n B" (o) from integrating over the semicircular path 


around the origin is 


7 [R v] enl + 


where Ro is the 


lim ou Han G(ké-yé) 1/2) (D.7) 


u>0 


R_(a) , where 


The remaining poles will contribute + | k 
1 


μυ] 


η 
R (a) = Tim (u-p,) S bun 6((42-02)1/2)7 (antr + 2) 


Uy Py k=- u 


where p, is the nth ordered pole of C(u) on the real axis. 
We may calculate the integral along the straight line portions 


of the path as 
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K(u) 2n[1 + T du 


620 k -u 
where 

Κυ) = 2 - 377 [L(u) + L(ue!”)] (D.8) 
and 

L(u) = $- [en 6 (ke -u?) '/4)) (D.9) 


The integral is interpreted as a principal value type, and it is 
denoted by the bar on the integral sign. The principal value of the 
logarithmic function is understood to have been used throughout. 


In summary, we have found that 





gn Bata) = K(u) 2n[1 + = ] du + HR =v) ipie 
δ»0 k^-u 


RM. ES α = α 
+5) Ría) - ) 2£m[1+ =]+)  2n[1 + —] 
e " nz] Pn m= | on 


and the value of G (o) can be written as 


R /2 00 00 
&(-/e0) «D? mn +21 0+2 
n=] Pn m=] m 


exp ju en[ ont ] - UA a Ϊ Cl == Pu 
s>0* E 
n 
. ο] X (D.9.1) 


n 
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If we perform a similar derivation in which we use B (a) 


instead of βία), as: the starting point, the result would be: 


gn B (a) = K(u) 2n[1 + 
620* k“-u 





] du * HR =v) 2n[1+ y] 


yn 
πμ 


(a) + 5 οι -Ἱ- Y anf] - E 
n=] Pn ΠῈ ] m 


Similar to (D.9.1), it follows that 


eoo 


] R /? - 
6 (a) » 80) (1 * D) ης η 


2 





n=] n me ] m 
ihy en ) - ὦν 
x exp[ = LE = ek * f K(u) 2n[1+ ]du 
$20 k -u 
j ? 


From (D.9.1) and (D.9.2), and knowing that Dum - G (a) , we 
recognize that the factorization is: 
R /2 © E 
6 (a) = 6-0) = 60) (1+8% n (1+ S97! 
n=] Pn 


x ia + E) exp | E anf tla- k qe din 
m=] eae 


oo 


t Ε K(u) 2n [ 1 * Siege, ΕΣ y lal] (D.10) 


$50" k^- a 
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To complete the theoretical portion of this appendix, we cite a 
theorem which provides for the expansion of a class of functions into 


infinite products [26,p.136]. 


Theorem: Let f(a) be a function with its only zeros at ἂ ἂν ἂν 
where lim a. is infinite and let f(a) be analytic for all a 
Then Fa) can be expressed 

f'(0 


a oo 
f(a) = f(0) e f 1 (1-Lje "7; 270 (D.11) 


If f(a) is both even and entire with its only zeros occurring at 
i - ΠΠ ` 
ta, ,tà,,fà4,..- and M tapi then f(a) can be written 


a f' (0 > SA. aa 
fla) = F(0) e * Te en.) 
EN = 


The exponential factors for the infinite products ensure convergence 
of the separate products. 


Consider 


a za 22 
μα νο ee ο... in 
n n 


Thus, 
za 
oo A co 22 
j Q- =) a je ) ed as noo 
n=] n ηΞ] η 7 


eo 2 
Because ) > converges as n+ , the infinite product also con- 
n=] n 
verges as n > 


co 9 
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We will now apply the foregoing factorization formulas to factor 


the expressions: 


(1,62) [K lb) = K lva) 77 


and 


[K (vb) [I (vb) - 1,(ra)]1 17 


Taking the top expression first, we define 


HD ia) = I (va) = 2.(a Vr- a) 
Dia) - x (b) - K (va) = i HHC) (b Vka) 


= MCA -a )] 


where use has been made of the well known relation between the modified 
and the ordinary Bessel functions [20,p.370]. 


It is readily seen from (D.12) that 


(1) e _ laa 3 E iaa 
f - J (ka) m (1-45 m qp (Ier—Je "" 
(a) a a) oi ( TT T^* I Πατ 


[z 
where a, = tilo, | ; lo, | = (De k* ; z, is the nth ordered zero 


of Jo) and τι - as noo , The value of ka must 


be restricted so that la | #0, therefore take ka << 1. 
Theorem C is used to factor eU (a). From the asymptotic 


value of the modified Bessel functions, we have 
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Uta. ~ ase se. [zi e -a a 


= a 2 cala] Bub'^ e/a " (318 


i a> © 


Thus ν- -5 and h=a. From (0.9), 


LU) » Scten e( C00) V2). Spon i Sen (ub) -uC (ua) 


o -b HET) (ub) +a HC (ua) 
"unm - uU (ua) 


From (D.7) 


R = Jim u le, = 0 L(u) has no other poles on 
2o OMNIS du 
the real axis, so Ra =0 


forall n. 


By (D.8) and the analytic continuation HZ) (zei) gim HET) (2) 


(1) (1) 
Di) i à. sir [m (ub) + aH (ua) 
HCT) (ub) E Ht ) (ua) 


-bH (2) (up) + at!) (ya) | 


HD (yyy =H) 


a 77 
Tl T 
A (ua) 


+ 2b[Y, (ub)J,(ua) - J, (ub)Y, (ua)] + 2a[Y (ua), (ub) - J,(ua)Y (ub)] 
.. AA IA IIA A O στ ο σσ 
[ 


J (W) - J,(ua) J" + LY (ub) - Y (ua)] 
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Factorization formula (D.10) gives 


gU (s) (Pito) - ἠ(])(κα) 1/3 εχρ[᾿Ὑ ρη[χία- 1 


KN (4) tapis Ja (D.13) 
620* k^- u 


X(o) * X (a) X (a) = {1 (va)[k (vb) - K (va)] } 


Thus 
x(a) 
X (a) = e ms (D.14) 
(1)4 " a 1αᾶ/ηπ 
G ^'(a) VJ. (ka) m (1* je 
ο pz To 
εχία) is introduced in the "plus" function and e~Xta) in the "minus" 


function to ensure algebraic behavior of each of these functions as 
jal +. Note also that y(-a) = -x(a) since Y (-a) = y (a) . We 
determine the value of y(a) from the asymptotic behavior of the 
denominator of equation (D.14). 


By (D.5) we know the asymptotic behavior of ο is 


2a 
k 








(1)4 1/4 ¡Gn 
ο. αν αἱ 5 as |α * 9 


in the upper half plane. The asymptotic behavior of the infinite 
product in the denominator of (D.14) can be found from the well known 


property of the gamma function [4,p.41], 
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b 
co r(— + J) 
- -Ca/a a 
"n mt a re a/an _ d 
n=] an +b NE S Om) 


Introduce Sterling's formula, 


T 
Ma) ve” a ER (29! 4 ...] as qo, 


larg a| « m 


Combining these two relations gives the useful result 


a b ] b 
co - — r(=+1) -[5 * =] 
lim I [! s ¿e lei eta 


Jom a 


[α]-ο ηΞ] 


x exp (0 - C - £n ΕΣ (p. 15) 


In (D.15) C is Euler's constant (0.5772 . . . ). Applying this re- 


sult we find the asymptotic behavior of the infinite product in (D.14) 


to be 
& 10d -ᾱ 
Tim Iu pues je Ta a een - c sogno + 47] 
ανν πι Παν ο ανω τὴ 
|arg al<n 


1, π 
where ο (na q) as noo, 


By collecting the above results, we determine that 











[α) - αἃ {1 - ς - gn|98| . 17 a 
1/2 


and that x (a) - x (-a) behaves as [al as |aļ >œ in the 


region of the a plane with 1 > -Κρ : 
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Turning to the other expression to be factored, 


K (vb) C1 (vb) - I (va)] 


we define 
2d (a) = 1, Gb) - 1,(va) - 4. (6 ψκξ-αὖ) - (a ia 
ka) = Kr) = i FAO) ita) 


Using (D.12), we have 


(2) m) 
Ν A Ber mm 
(a) LJ, (kb) - J 9 (ka) ] m TT T^ 
2: 1a atb) 
x I (1 + α MT a 
m=] Τα" 


where an tila, | , 


| 2m 2 
[74 E (=) - k* and zo is the mth ordered zero 


of Jo) - J ($ z) 


The value of kb is restricted to kb «« ] , so that || 10. 

The asymptotic behavior of the zeros of [0 (2) - J (E z) are 
determinable by considering the asymptotic behavior of the Besse] 
function for large arguments and considering that the difference be- 
tween the values of a and b does not affect the magnitude. Differ- 


ences between the values of a and b , however, cannot be neglected 
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in phase factors. Through these considerations and trigonometry, we 


find that the most rapidly varying roots are those given by 


2 + (m- y 5m as moo 3; 
nc IBS. 


Since 


aa σας erblal as ato 


then v= - 5 and h=b 


(1) 
(2) mud semen UN) » sen (ub) 
L (u) AT [anti 7) Ho (ub )] = FUE 


a 
nn 
N 
u 
— 
E 
u 
ajo 
ze 
σ 


rn): 2 ] ] 
D mub J^ (ub) + Cub) 


CE [5 i HET (xp) 31/2 exp| 2. sa [E18 ] - e 


+ f Ku) n + eL du | (D.17) 
-0t Y 


6 k`- 


Collecting the results, we get 
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Yla) = Yla) Yla) = {KIEL Er) = Ira)” 


where 


e $9) 


Y” (a) = c Ta 
A - Jai ind iz 
ο ο me] ilo | 


(D.18) 


By (D.5), we know that 


- «Ντο 
emo ela. d e. eges 


in the upper half plane and from the asymptotic behavior of infinite 
products (D.15), we know 


atb 








a ) 1 
e Ta * ToT ν e MIO Se 1-C 
| arg a|«m 
alatb TW 
- an | Bh e i53) 
2 
DEZ LEE Zoe 
Therefore 


in , 2b 20 
tu =p MIR ] (D.19) 


αί atb) 
T 


as lal >œ in T> -k . 


dla) = a ath cs 








ΠΝ ΠΠ 1ο). ]ε]| s 
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Appendix E - Radiation from the Infinite Cylindrical Antenna 


In this appendix we shall determine the radiation from a hollow 
infinite cylindrical antenna. An integral equation for the current on 
Such an antenna will be formulated and solved. The fields will then 
be determined from the current. The integral equation considered in 
this appendix was first formulated and solved by Hallén [35,36]. He 
was interested in the reflected current on a cylindrical antenna, when 
the incident current wave comes from a source that appears to be an 
infinite distance from the end. Although the integral equation is the 
same as that found by Hallén, the application used in this appendix is 
novel. 

Consider Figure E-1, which is an illustration of the actual 
geometry and its equivalent representation, which is amenable to 
analysis. Since the antenna is made of perfectly conducting material, 


its surface is an equipotential and 


ESSE) = -v(t) 6(z) when p = a, ~<2< 
In the —— domain this becomes 

E (Tw) = -v(w) 6(z) 
In view of the rotational symmetry, it can easily be shown that 


2 > 
A_(r 
E (T,u) = *iu[A (V) + 2] O) 


Y is the vector from the origin to a point on the surface of the 
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P(a,$,z) 


$ 
Lo 





Figure E-1. Infinite cylindrical antenna and equivalent repre- 
sentation using gap generator as source 
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antenna and A (T) is the well known vector magnetic potential, which 
satisfies the relation B=VxA 
The symmetry of the structure also requires that A, (T) = A (z); 
thus the axial component of the vector magnetic potential satisfies the 
differential equation 
Jen (z) 


Z 


SL 
s + Ka (z) = + E v(u) δία) (E.1) 
ðZ 
We also have from symmetry that A (z) = Α.(-Ζ). 
The solution to (E.1) obtained in a straightforward fashion is 
ik|z| 
_ f kv(u)e! 
A(z) = C cos kz + C sin kz + = - 
From the symmetry condition that A(z) = A(-z) and the requirement 


for only outgoing waves on the infinite structure, C, = C, - 0 and 


IK 
A(z) = kv(w)e KIzl_ (E.2) 


2ω 


The vector potential will now be determined by an alternative fashion, 
whereby use is made of the free space Green's function and the surface 


current. The general formula is: 


> +> 
E" ik| r-r' | 
Kcr) = 2 co cM. qe 
E [ΠΠ 


with 


Ir - 7 |= /(z-z')% Sena en 200 'cos(-¢') 
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In the case of the infinite cylindrical antenna with the obser- 


vation point on the cylindrical surface, this general formula becomes 





co 2n 
AA 
A (z) E | Hz) | E d$'dz (E.3) 
00 0 


whe re 


r= J(z-z')'* 4a' sin (4) 


and 6 is taken as zero. Since the vector potential is independent 
of  , the particular choice of is immaterial. Since (E.2) and 
(E.3) are both equivalent expressions for the vector potential, they 


may be set equal to each other. The results are 





oo 2m . 
ikr ik zul 
00 0 


27 MW 
e P 
= m (= v(w) eikizi (E.4) 
Ho 


where -o«z«o 

The equation (E.4) is valid for all values of z , since it is 
based on the mathematical model of the actual antenna structure. 
However, in the actual antenna we observe that the current on the 
feedline after t - 0 , the time when the incident excitation reaches 
the z= 0 plane, is a TEM mode wave. The current wave can be ex- 


pressed as: 


I(z) = 1(0) e IK i z<0 
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I(0) is the magnitude of the current at the plane z=0 . In view 
of the above, it is justified to represent the antenna feedline system 


with the following set of equations: 





€ 3 
oo 2T m 2T hae vlwJe KZ u τη 
στ | dome.) | e do! = 0 

0 


--ΟΟ 


g(z) | 2 <0 


le = ο) αρ πνς E.o È) 


We restrict the analysis to the consideration of TEM mode feedline cur- 
rent only. The function g(z) for z «0 is unknown. The solution to the 
integral equation (E.5) is the current on the infinite cylindrical 

antenna when there is an outward (from the origin) travelling current 

wave on the feedline. 


The current is now expressed in tems of its spatial] Fourier 


inverse 
Me) =a | eon ] (8 da 
ul 
I(a) = I, (a) + I (a) 
with R 
I (a) = Lf I(z') οἱ ασ dz ' 
ven 6 
0 — 
(a) = * I(z') e!97 gaz 
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-K5z ΚοΖ 
Note that |I(z')| «Ae s z>o wind |I(z')| <Be as 


zZz > -~ where Ko is the imaginary part of the propagation constant. 
I(a) is obviously analytic in the strip -k, ES Κο of the 
a=o+ it plane. I (a) is analytic in the half plane t > -k, and 
I (a) is analytic in the half plane rt < SE 

From (E.5) it follows that 


0 
I (a) = κα | I(0) LS. et laz d= 2 j πο) — 
ver Yon (a-k) 


ου τς Κρ . Additionally, 


G (a) - m m 
ven) |W 
En 


ο 
[2πε 
ρα e 10) V u) - 
= 1 ues hae for T > Κρ 


Upon putting these results back into (E.5), we obtain 


v(w) οἱ Κ7 οἱ αξ d 





Z 





oo 00 2T $ eo 
ji ikr Pe 
| dz' —— [ο I (o) du | — — | e^197 G(u) da 
M von S 0 ca 


or equivalently 


oo oo 


ἐπ . 
ter ο ο as 
an | a | ] us 7s (o > | r P 





oo 


= | eae (¡> B + G (a)] do (E.6) 
ο 


00 
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We digress at this point to find an equivalent form for 
oT 


i 
| — de‘ . It is well known that the free space scalar Green func- 
0 


tion in spherical coordinates satisfies (vés k*) Glee) = cigar); 
and ΤΝ 
ik[r-r'| 
- >, EINE 
G(r,r') = E - 


Ir-r'| 


In cylindrical coordinates the Green function satisfies 


UR A AI 


p 
and 


G(o,0',0-$',z-z') - —, } ¿im(9-6*) 


Anm 
x | ο σας’) I (ve_) K(ve,) da 


In the above equations y= C Oz and o, is the 
> 


be of p and p' . These two Green's functions are equiva- 


lent three-dimensional representations. Therefore, 


A de 
E rue . n YO gos 
and 
IT > > oo 
imer. | : ' 
e E -ia(z-z') 
| aa dp' = 2 [e (ro) Ko(vo,) do 
0 E -00 
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We now evaluate the above equation at po = 9' = a and sub- 


sertute Tt into (E.5). 


1 f pl -iaz'rp,* Pike f -iß(z-z') 
] dz Ji ο” α΄ [11 (ο) ο da | e P071) p Ga)k (ya) de 
a iv(w) , G”(a)] da (Επ) 
"y (atk) 


Interchanging the order of the integrations on the left hand 


side of the above equation, and recognizing that 


co 


| ertlaß)z" aj. - 2ηδία-β) 


«OD 


equation (E.7) transforms into: 


-182-,+ 11(0 - 
2 | da ] e782r7*(q) m 3 toes va) Ste) d 
= | e 2 [ To πο τς) G (a)] da 


Integration over the ß variable produces the result: 


2 | πος ([ (a) - mes [I (va) Koya) J 
πία-Κ 


TE x 
T (w) G (a) 2 
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Let I (ya) K bra) be factored into the product L" (a) L (a), 
where E) is analytic in the half plane qt > -k, and L (a) is 
analytic in the half plane rt < Κο ἢ 

4 


After substituting the factored form into (E.8) we derive the 


following result: 


2 | anto - „pe: haw 
«DO a 


3 m [i e VW) | G(o)] » 0 
L' (a) "o (atk) 


The integrand is regular in the strip k, er Ko of the plane. 


Since it must vanish, it gives the result 


LU. VO) 0) ue : v(w) 
v2n(a-k) οι (atk) L7 (a) 


Multiplication of both sides by a k* yields the integral or entire 


function 





L*(a) I*(a)(o°-k*) - HAMM o = Gla) (42.42) 
T 


eL (a) 
+i ¡o v(w)(a-k) _ Ρ(α) (E.9) 
Ho L (a) 


The left hand side is analytic in the upper half plane (1 > -Κο}, 
the right hand side is analytic in the lower half plane (t < Ko) and 


each side is the analytic continuation of the other. We divert our 
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attention from (E.9) for the moment to consider Io); L (a) and 


their asymptotic behavior as a +œ in their respective half planes 
of regularity. 


Earlier, we defined the functions 
+ E 
L (a) t (a) » Ij(va) K (va) 


By the methods of Appendix D, this factorization can be carried out 


with the results 





ES a) vJ (ka) i (re ™ 
συ ο ---ω. 
RO 


g^ 3 a) = [> int D (ka)] V2 exp| 12 E el ] - 


+ KG) (uy gn[ I ‚Ale η du | 
620^ k`- u 
ae, „2 m 
E RA y P neos (ua) * Y^ (ua) 
0 0 
Dp 
via) = = Q-C- 1π| S1 + an[ Slt is] 
2m2_ 12 
Ιδίως = - k where z_ is the mth ordered zero of 


y, (2) 





-177- 


Lla) « ee as Jal »® in half plane r > -k . 


wO j 7° 


as lal » in half plane τς Κρ. 

Before returning to (E.9), we note from the edge conditions 
[21] that the current goes to zero at the terminus of the hollow con- 
ductor. As I' (o) and G (a) represent transforms of current we can 
deduce the asymptotic behavior of these functions for a by 
examining the behavior of the respective currents as we approach the 
ends of the antenna. By considering the antenna length to be finite 
and changing the variable to place the origin at the terminal ends, 
it is easily shown that I (o) and G (a) vanish as la] >> in 
their respective half planes of analyticity. 

„Application of the asymptotic behavior of the various func- 


tions to (E.9) shows that 


< lal as a>% , 





Mo ὑπ = i μία 





TSS -ks 
and 


κ. o v(w)(a-k 
e L (a) 








ι΄ as ato , 





T < Κο 


From the extended form of Liouville's theorem, Pla) is a polynomial 
of degree less than or equal to one [4, p.38]. We may therefore write 


P(a) = C,+ C (a-k), where C. and C, are constants to be determined. 


1 2 
By setting a= k in the left side of equation (E.9), we get 
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ο πο 


τ 
Additionally, 
2 K $ iI(0) L'(a) (o*k) 
— + (0 = La) Tla)lark) + ee 


Yen (a-k) 


By letting a >+* in the above expression, we find C, = 0 . Thus 


Pla) is a constant. 
From (E.9) and the value of P{a) just determined, we have 


tay - H0)... i21(0)k L'(k) 
ven (a-k) a. jL (a) 


Recall that I(a) = Ee + I (a), so the transfom of the current is 
, + 
(a) - . 121(0)k τ (6) 
Jana) (a) 
_ . i21(0)k L’(k) L”(o) (£.10) 
ven (a -k 


JI, (va)K (a) 


Inverse Fourier transformation of the above yields 


oo 


(a) - - Ἠ(ὼκ MOS f εις i (t A 


(a -k )I S (va)K (va) 


oD 


It is also clear that if we close the contour in the upper half plane, 


the integration produces I(z) = (oje ikz, z <0 , as required. 
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The vector magnetic potential can now be found from the general 


formula, 
οο 2m AF >, 
O πλω eikir-r| a. 
Α = — ————— —— d$' dz 
Z ar T >>, 
-00 0 [r-r | 


BR Uk) ik FF 
A | | | dz! @ 
2 4? 27 


-0 -o0 | r-r' | 


J al az' L" (a) m 


(a-k) J (ya)K (ya) 


oo 


iu I (0)k L'(k) 


ONE — E | ο κος, ) 
T 


=00 -Ὁ 


LN paw α- 
Κ (ο Je -k )] | e L a) da 
L (af-kê)1 (va)K la) 


iu L(0)k L*(k) Í do ae Rk (ο Je? -k9) L laje 182 


- ---η----- Do '.'"-".. = | ---------- 
2T κ. (a^ - k ) I (ya) Kota) 


x 6(a-8) 


eoo 


| seek (yo) | La) 
e i). a 
(a =k") K (ya) 


iu 1(0)k L'(k) 


= 
QT 


=O 


This equation is valid for p>a,2z2>0. Note that the magnetic 


flux density is given by 





therefore the magnetic intensity is 


il(Q)k t*(k) % e7taz Kılyp) L'(a) da 
A, (052,0) = | NN (E.12) 


Equation (E.12) is exact for an infinite length antenna. If we assume 
that it also holds for an extremely long antenna and employ the 
methods of steepest descent aS was done in Chapter 4, the resulting 


equation for the magnetic intensity is given by 


_ I(0)e'kr J (ka sin 8) f k) 
H¿ (0,230) > A πο EP (EMS) 


Earlier we assumed that I(0) was the incident current which is 
reflected from the discontinuity at the Ζ Ξ 0 plane [8]. The re- 
flected current on the feedline from pulse excitation is known through 
experiment to approximate a reflected incident current pulse. From 
the edge conditions for this structure it is clear that the current 

on the outer conductor vanishes at z 720 . If the current on the 
outer and inner conductors are to be equal in magnitude and opposite 
in direction, then the current on the inner conductor must also vanish 
at 2720 , the point of discontinuity. We therefore conclude that a 
pulse of current is launched onto the antenna and an oppositely 
traveling pulse is launched back upon the center conductor of the feed 
line. The magnitude of this current is assumed to equal the magnitude 


of the incident pulse. 
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If the incident field is given by equation (2.7) of Chapter 2, 


then 
I(0) - -2π fo uy 
-m αν u 
o m3 


The negative sign is required since the reflected current travels in a 
direction opposite to that of the incident current. [Incorporation of 


the above into (E.13) gives the final form of the magnetic intensity, 


ikr J „(ka sin 6) r 
Al E auc A 
$ L (k cos 0) 
rf 


+ 
lim SO = | 


ka* 0 L (k cos 0) 


then (E.14) obviously reduces to 


a... fo M elkr J (ka sin 8) T 
oet Ho nÈ r” sin 6 
a 


for an antenna with kh >> 1 and ka << 1. 
To show the validity of this limit we rewrite it from the ex- 


pressions for L*(k) and L (k cos 0) , 


ika 
00 k ET 
T E aB pe eG) RU k cos 0) 
"NL A 10 BOà^31 .Á . , 
ka» 0 A Bi ika cos 0)/ (mr) (AS a SUO 


m=] 





= > > => 
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exp KG (uy £n[ 1 + E ]du 
8-0" k -u 
= Jim 
ka +0 E 
exp Ϊ Κέι) 2n[ 1 + BO y 
TM ow c 


= lim exp Ϊ K*(u) gn[1 + a 39 du 
>) 


kan SUM cos 0 * Jj] - (u/k)* 


This final limit was investigated in Chapter 5 where it was shown to be 


uni ty. 
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